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ABSTRACT 

We use a novel method to predict the contribution of normal star-forming galaxies, merger- 
induced bursts, and obscured AGN, to IR luminosity functions (LFs) and global SFR densities. 
We use empirical halo occupation constraints to populate halos with galaxies and determine 
the distribution of normal and merging galaxies. Each system can then be associated with 
high-resolution hydrodynamic simulations. We predict the distribution of observed luminosi- 
ties and SFRs, from different galaxy classes, as a function of redshift from z = — 6. We 
provide fitting functions for the predicted LFs, quantify the uncertainties, and compare with 
observations. At all redshifts, 'normal' galaxies dominate the LF at moderate luminosities 
~ L* (the 'knee'). Merger-induced bursts increasingly dominate at L S> L*; at the most ex- 
treme luminosities, AGN are important. However, all populations increase in luminosity at 
higher redshifts, owing to increasing gas fractions. Thus the 'transition luminosity' between 
normal and merger-dominated sources increases from the LIRG-ULIRG threshold at z ~ to 
bright Hyper-LIRG thresholds at z ~ 2. The transition to dominance by obscured AGN evolves 
similarly, at factor of several higher Lir. At all redshifts, non-merging systems dominate the 
total luminosity /SFR density, with merger-induced bursts constituting ~ 5 — 10% and AGN 
~ 1 —5%. Bursts contribute little to scatter in the SFR-stellar mass relation. In fact, many 
systems identified as 'ongoing' mergers will be forming stars in their 'normal' (non-burst) 
mode. Counting this as 'merger-induced' star formation leads to a stronger apparent redshift 
evolution in the contribution of mergers to the SFR density. We quantify how the evolution 
in LFs depends on evolution in galaxy gas fractions, merger rates, and possible evolution in 
the Schmidt-Kennicutt relation. We discuss areas where more detailed study, with full radia- 
tive transfer treatment of complex three-dimensional clumpy geometries in mixed AGN-star 
forming systems, is necessary. 

Key words: galaxies: formation — galaxies: evolution — galaxies: active — star formation: 
general — cosmology: theory 



1 INTRODUCTION 

Understanding the global star-formation history of the Universe re- 
mains an important unresolved goal in cosmology. In recent years, 
observations of the properties of galaxies in the infrared, at red- 
shifts z = — 3, have begun to shed light on this history, but have 
also revealed a number of intriguing questions. 

Of particular interest are the roles of mergers and AGN in 
driving star formation and/or the infrared luminosities of massive 



* E-mail:phopkins @ astro.berkeley.edu 



systems. A wide range of observed phenomena support the view 
that gas-rich 1 mergers are important to galaxy evolution; but it 
is less clear what their role is in the global star formation pro- 
cess and buildup of stellar mass in the Universe. In the local Uni- 
verse, the population of star-forming galaxies appears to transition 
from "quiescent" (non-disturbed)" disks - which dominate the to- 

1 By "gas," we refer specifically to cold, star-forming gas in galaxy disks, 
as opposed to hot, virialized gas. 

2 In this paper, we use the term "quiescent" to refer to star-forming systems 
that are not strongly disturbed in e.g. major mergers and are forming stars in 
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tal star formation rate/IR luminosity density - at the luminous in- 
frared galaxy (LIRG) threshold 1O U L (M» ~ lO-2OM yr _1 ) 
to clearly merging, violently disturbed systems at a few times this 
luminosity. The most intense starbursts at z = 0, ultraluminous in- 
frared galaxies (ULIR Gs; Lir > 10 Lf^), are invariably associ- 
ated with mergers (e.g. | joseph & Wrigh t 1985; S anders & Mirabell 
1 19961 : lEvans et alj|2009b~ with dense gas in their centers provid- 
ing material to feed black hole (BH) growth and to boost the con- 
centration and central pha se space density of merging spirals to 
match those of ellipticals dHemquist et alj [l 993; Ro bertson et al.l 
l2006h . Various studies have shown that the mass involved in these 
starburst events is critical to explain the relations between spirals, 
mergers, and ellipticals, and has a dramatic impact on the proper- 
ties of merger remnants (e.g ., Lake & D ressle dl 19861; Dovon et al. 



1994l:IShier & Fisched[l99lllJames et alll999HGenzel et alj|200ll; 



Tacconi et al 



2002; Dasyra et al. 2006, 2007: Rothberg & Josep 



2004l2006l : lHopkins et alj|2009aft eb~ 



At high redshifts, the role of mergers is less clear. It is clear 
that LIRGs and ULIRGs increase in relative importance with red- 
shift, with LIRGs dominating the star formation rate/IR lumi- 
nosity densities at z ~ 1 and ULIRGs dominating a t z ~ 2 (e.g. 
Le Floc'h et all 120051; [Perez-Gonzalez et alj 120051 ; ICaputi et al.1 
2007l : lMagnelli et alj2009b . This, together with the fact that merger 



rates are expected and observed to increase with r edshift (by a fac- 
tor ~ 10 from z = — 2; see e.g. iHopkins et al. I l2009il and refer- 
ences therein) has led to speculation that the merger rate evolution 
may in fact drive the observed evolution in the cosmic SFR den- 
sity, which rises rapidly from z ~ — 2 and then turns over, de- 
clining more slowly (e.g. lHopkins & Beacom|[2006l and references 
therein). 

However, many LIRGs at z ~ 1, and potentially ULIRGs at 
z ~ 2, appear to be "normal" galaxies, without dramatic morpho- 
logical disturbances associated with t he local starbur s t population 
or large apparent AGN contributions dYan et al] [2007; Saii na et all 
l2007l : iDevetal]|2008! ; iMelbourne etai]|2008l ; iDasvra et alj|2008b . 
At the same time, even more luminous systems appear, includ- 
ing large populations of Hyper-LIRG (Hy LIRG; Lir > 10 B L n ) 
and bright sub-millimeter galaxies (e.g. | Chapman et alj 1 20051 ; 
I Younger et ai]|2007l . l2009d : ICasev et alj|2009b . These systems ex- 
hibit many of the traits more commonly associated with merger- 
driven starbursts, including morphological disturbances, and may 
be linked to the emergence of massive, quenched (non star- 
forming), compact e l lipticals at times as early as z ~ 2 — 4 



(Papovich et al. 2005; Younger et al. 2008bk iTacconi et al 
ISchinnerer et alJl2008l : Tchapman et alj|2009t ITacconi et al 



120061 : 



2008 



But reproducing their abundance and luminositie s remains a 
challenge for current models of galaxy formation |Baugh et al 
2005l:ISwinbank et aLbOollNaravanan et alj2009b] ; l Younger et~ 



2009al) . 

In a related vein, observations of a tight correlation between 
the masses of super-massive BHs and their host spheroid properties 



(Gebhardt et ; 



super-m 
ail hood: 



sp here 

Ferrarese & Merritt 2000l : Magorria n et al.1 



ll998l : lNovak et alj2006l ; lHopkins et al.l2007bl) suggest a tight 
pling between BH growth and star formation, perhaps in particular 
to the mergers believed to drive the formation of the most massive 
bulges. Cons idering the en ergy output required to form the BH pop- 
ulation (e.g. ISoltarjl 1982h . or the obser ved bolometric quasar en- 
ergy density as a function of redshift (see lHopkins et al ]2007d and 



similar fashion to most "normal" disks. We do not mean non-star forming 
systems, as the term is used in some literature. 



references therein), it is clear that the bolometric output of quasars 
and AGN is at least roughly comparable to the total infrared lu- 
minosity density of the Universe at most redshifts (z ~ 1 — 3) — 
although the measurements above suggest it is still a factor ~ 2 — 3 
lower. Some recent observations have suggested that the population 
of very luminous, highly obscured (Compton-thick) quasars may 
be considerably larger than previous estimates, in which case the 
heavily-obscured AGN population could represent a large fraction 
of the total IR luminos i ty density at high r edshifts dHickox et al.1 
l2007l : lDaddi et al.l2007l : lTreister et al]2009b . This would have dra- 
matic implications not just for BH populations and e.g. the implied 
radiative efficiencies of BH accretion, but also for the implied total 
star formation rate density. Some apparent discrepancies between 
e.g. the total mass density observed in old stars and the implied 
star formation rate density have been cited as possible evidence 
of a time-dependent stellar initial mass function (IMF); but a rising 
contri bution from obscured AGN at high reds hifts could mimic this 
effect dHopkins & Beacorrtl2006l : lDave1l2008| ). 

In particular, there are long-standing questions of what pow- 
ers the most luminous infrared sources, for example, ULIRGs and 
sub-millimeter galaxi es. This debate extends to the discovery of 
these objects (see e g.lSoifer et al J 1 9841.1 1 9871; IScoville et al.ll98rj ; 



ISargent et al.ll987l : ISanders et al.ri988bl : ISolomon et al] 1997b . and 
has persisted despite the addition of millimeter spectroscopy and 
observations in a large number of independent wavebands (for a 
review of the debate, see both ISandersj[l999l : Ijosephl 1 1999b . Al- 
though some evi dence suggests tha t they are prima r ily powered 
by star formation dF^rah et alJl2003l:lLutz et all 19981: ISaiina et all 
l2007l : |Pope et al]2008j |bt lWatabe et al.ll2009l : lNardini et al]2009b . 
the constraints and correlations typically invoked have inherent fac- 
tor ~ 2 uncertainties, and thus could easily accommodate compa- 
rable power input from star formation and AGN. Moreover, a suf- 
ficiently obscured AGN, in a medium with the right optical depth 
properties, is indistinguishable from star formation by the usual in- 
dicators (e.g. PAH strengths, emission region sizes, or any other 
infrared spectral or morphological criteria). Hence even at z = 0, 
debate surrounds the power source of many bright infrared systems, 
and there exist a number of examples of systems classified as "star 
formation dominated" by all of these metrics that later revealed 
Compton-thick AGN whose longer-wavelength e mission has been 
fully re-processed, even into "cool" dust (see e.g. lAlexander et all 
l2008l and references therein). In a bolometric sense, the most lu- 
minous galaxies observed (with L > 10 14 Lq or ^> 10 47 ergs _1 ) 
are the most luminous quasars; although the contribution of these 
systems to the infrared remains highly uncertain. This may be im- 
portant for resolving the theoretical difficulties in modeling these 
bright systems. 



ing these processes in an a priori manner (see e. 


g. iBaugh et al. 


2005; 


Hopkins et al]2008d; Narayanan et al]2009b 


lYounger et al. 


2009a 


). However, two basic limitations remain. In direct cosmo- 



logical hydrodynamic simulations, as well as semi-analytic models 
of galaxy formation, it is well known that it remains challenging 
to accurately reproduce global quantities such as the galaxy mass 
function and the distribution of sizes, gas fractions, and hence star 
formation rates, especially the distributions of star-forming gas and 
their rela tions to whether or not g alaxies are "quench e d" (recently, 



see e . g. IWeinmann et al.| 
120081 : 



Fontanot et al 



2009) 



2006; Mail er et alj boOfj ; iKimm et al.1 



This makes it difficult to determine 
whether discrepancies between such models and the observations 
owe to their treatment of star formation, or to discrepancies in 
these quantities. Moreover, it is difficult to disentangle the effects 
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of these different properties on the distribution of star formation 
rates. In addition, for merger-induced starburst and AGN activity, 
although it may be possible to roughly estimate some global quan- 
tities (e.g. the total mass involved in a starburst) from simple ana- 
lytic motivations or low-resolution cosmological simulations (sev- 
eral ~kpc typical), it is not straightforward to estimate the chaotic, 
time-dependent behavior of full lightcurves needed to estimate the 
distribution of time spent at different, rapidly varying luminosities 
without high-resolution simulations of individual systems. Since 
the number density of the most bright systems is exponentially de- 
clining, fluctuations and features in the starburst/ AGN fueling his- 
tory on small time and spatial scales (At ~ 10 7 yr, R < lOOpc) can 
be critical for correct estimates of their contributions to bright pop- 
ulations. 

In this paper, we present theoretical predictions for the distri- 
bution of galaxy star formation rates and infrared luminosities, as 
a function of galaxy mass and redshift, using a novel methodology 
that can circumvent some of these obstacles. We combine a halo- 
occupation based approach, in which we take galaxy properties as 
fixed from observations at each epoch, and then apply rules for the 
distribution of star formation rates/infrared luminosities in "quies- 
cent" systems, merger-induced starbursts, and obscured AGN, cal- 
culated from a large suite of high-resolution hydrodynamic simula- 
tions of individual galaxies and galaxy mergers. We use this to inde- 
pendently estimate the contributions of "normal" galaxies, mergers, 
and AGN to the luminosity functions. The comparisons we make 
are approximate - we do not include full time-dependent radiative 
transfer in simulations (the subject of future work, in progress), and 
so focus on integral quantities such as the total IR luminosity and 
SFR distributions, that are less sensitive to issues of e.g. the exact 
dust distribution, temperature, and other properties. We also explic- 
itly separate the contributions of AGN and star formation, but stress 
that, in real systems where the two are comparable, their additive 
effects are non-linear, and will require further study. We show how 
adding or removing components of the model taken from obser- 
vations such as e.g. the distribution of galaxy sizes and gas frac- 
tions affects these consequences. We compare to observations of 
all quantities, where available, and find reasonable agreement but 
with some interesting apparent discrepancies at high redshifts. We 
also show how these populations relate to the scatter in star for- 
mation rates at fixed galaxy masses, and in a global sense to the 
total star formation rate density, the star formation rate density in 
mergers, and the fraction of the inferred star formation rate density 
which might really be driven by obscured AGN activity. Readers 
interested primarily in the comparison of predictions and observa- 
tions may wish to skip directly to §[3] 

Thro ughout, w e adopt a JIm = 0.3, Qa — 0.7, h = 0.7 cosmol- 
ogy and a lChabrierl d2003l) stellar IMF (discussed further below), 
but these choices do not affect our conclusions. 



2 THE MODEL 

The model used here is a slightly modified version of one that 
has been discussed e xtensively in a series of papers (most recently 
iHopkins et al. We summarize the salient properties here. 

FigureQ]provides a simple outline of the model, on which we elab- 
orate below. 



2.1 Halo Occupation Constraints: The Initial Galaxy 
Population from Observations 

At a given redshift, we use the halo occupation distribution (hence- 
forth HOD) formalism to construct a mock sample of galaxies. 
Specifically, we begin with the observed galaxy stellar mass func- 
tion (MF), which we take as given. Since we are interested specifi- 
cally in star-forming galaxies, we adopt just the galaxy stellar mass 
function of star forming or "blue" galaxies where available (here 
at all z < 2); although there may be some trace star formation in 
red galaxies, assuming typical values yields a negligible contribu- 
tion to the bright far-IR and massively star-forming populations. 
At redshifts z > 2, type- separated MFs are no longer available, so 
we simply adopt the total galaxy MF (i.e. assume all systems are 
star-forming); however, the fraction of massive galaxies that are 
"quenched" and red has become sufficiently low by z = 2 (and is 
rapidly falling) that it makes little difference (e.g. adopting the up- 
per limit - that the red fraction at all masses at z > 2 is equal to that 
at z — 2 - makes no difference to our predictions). 

The uncertainties in the galaxy abundance are one of the dom- 
inant uncertainties in the model, especially at high redshifts. We 
therefore consider different mass function fits, to represent the pos- 
sible range. At z = 0, the uncertainties are relatively small; we 
adopt the mass function of star forming galaxies from iBell et al.l 
j2003l) . Fromz = 0-2, we con sider the MFs of s tar-fo rming galax- 
ies from lArnouts et alj J2007h and lllbert etal] d2009h . The range 
between the two is representative of the uncertainties and scat- 
ter in a number of other calcul ations, which cove r different por- 
tions of this dynamic ranee (e 
20061; iFranceschini et aimOO 



Bundv et alj 



Borch et aU 



20051; IPannella etal 



20061 ; iFontana et al 



200d ; lBrown et alj2007l) . At z > 2, we adopt as bracketing the rele- 



vant d yna mic range the mass fu nctions from lPerez-Gonzalez et alj 
J2008I) and Fontana et al. (2006); again, other determinations (e.g. 
iMarchesini et alj2009l ; lKaiisawa et alj2009l) lie within this range. 

Given each galaxy and its stellar mass, we assign it other prop- 
erties in accord with observations. First, a gas mass. It is well- 
established that, at fixed stellar mass, galaxy gas fractions are 
higher at high redshifts (see compiled references below). More- 
over, the trend of galaxy gas fractions as a function of stel- 
lar mass, and their scatter, have been quantified (both directly 
and indirectly) at a range of redshifts from z — to z = 3. 
We have compiled observations from the available sources, span- 
ning this redshift range and a stellar mass range from M« ~ 
10 10 — 10 i2 Mq (more than sufficient dyn amic range for t he pre - 
dictions of i nteres t here), specific al ly fro m[Bell & de Jong (2000); 
McGaughl d2005h : Icalura et alj <|2008|); IShaplev et al.l (|2005 b ; 



Erb et alj d2006l); IPuech et alj J2008h; uMannucci et alj h()0% : 
Cresci et al.N2009h;lForster Schreiber et alJj2009l);lErbl J2008h. We 



presen t t hese observati ons in a number of papers ( IHopkins et alj 
l2008dl lbl 120073, 12009bl l3) and show that the z = zero-point and 
evolution with redshift can be well-fitted by the simple functions 

/-(M>|t = 0)=Ji, W l + (JM>/ 1 10Mayu 

/..(M. |z) = /o{l- r(z) [1 - fo /2 ir 2/3 , (1) 

where r(z) is the fractional lookback time to a redshift z (= at 
z — and = 1 at z — » oo). The former functional f orm is mo- 
tivated by cosmological hydrodynamic simulations ("K eres et alj 
I2005L I2009T) . and the latter by the sca lings of simple closed-box 
models that obey the iKennicut t Nl998h relation at all times. " The 

- This function is presented in IHopkins et aU l2009ch . Equation (2); we 
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Input/Initial Conditions: 



Adopt Observed Galaxy 
Stellar Mass Function: 
n(M 5tars 



Assign each galaxy a size and 
gas mass from observed 
correlations: 

stars J 



r 



Use observed abundances 
& clustering to populate 
HOD (place each galaxy in 
a halo). 



Knowing the halo 
evolution, this gives the 
>alaxies that have recently 
had a merger: 

f m erger( M stars, ) 



Add Models: 



Determine starburst mass Mburst & 
peak BH luminosity L pe ak from high- 
resolution merger simulations, as 
function of ( /t, f £ls , R e , M stars ) 



Construct Luminosity 
Functions: 




c 



Adopt observed 
obscured fraction and 
template IR spectra, as 
function of Lbh 




Star 
Formation I 




J 



Assign "steady-state" SFR from observed 
Kennicutt-Schmidt relation: 
2sFR~2 gas L4 ~(M gas /Re 2 ) 1 - 4 



Figure 1. Summary of our model methodology, described in detail in §[2] We begin with a halo occupation model: at a given redshift, the galaxy stellar mass 
function, and distribution of galaxy sizes and gas fractions are taken from observations. Placing galaxies in halos from dark matter simulations according to their 
observed abundance and clustering, and evolving these forward in time some short time, we obtain the merger population. Non-merging galaxies are assigned 
a star formation rate based on their size and gas mass, according to the observed Kennicutt-Schmidt relation; for mergers, the lightcurves of merger-induced 
starburst and AGN activity are taken from fits to high-resolution galaxy merger simulations (including star formation based on a local Kennicutt-Schmidt 
law, gas cooling, and feedback from accretion and star formation), as a function of salient galaxy properties. Star formation rates are converted to infrared 
luminosities with a simple empirical proportionality; AGN bolometric luminosities are corrected to far-IR luminosities based on observed obscured fractions. 



important quantity is not the precise scaling, but rather, the fact 
that it provides a convenient interpolation formula between the ob- 
servations above. We assume a constant, intrinsic 0.25 dex scatter 
about these gas masses at each stellar mass, also in agreement with 
the observations above. 

Next, we assign galaxy spatial sizes, again from observa- 
tions. At z = 0, the distributions of disk sizes are well-measured; 
at higher redshifts, there is some uncertainty, but observations 
are converging on the conclusion that the star-forming population 
evolves relatively mildly in size with redshift (whereas the non- 
star-forming population evolves more rapidly). From z — — 2, 
we find that the compilation of observation al results on t he evo - 
lution of the disk s ize-m as s relation from iTruiillo et alj (120041); 
iRavindranath et al.l J2004I); iFerguson et alj 1 20041) : iBarden et all 
d2005l) : iToft et alj J2007I): lAkiyama et all J2008h . and the theoret- 
ical models in lSomerville et alj J2008h . can be simply represented 



note that in the text there is a typo, and the equation is written with a +t(z). 
The correct form (above), with — t(z) was, however, used for the calcula- 
tions therein. 



as relatively weak power-law evolution in disk size at fixed mass 



R e (M* \z = 0)^5.28 kpc( 



10 10 M 



© 



R e (M*\z)=R e (M* |z = 0)(l+zr ' 6 (2) 

where the z = relation is taken from lShen etal] J2003T) (appro- 
priately normalized for our adopted cosmology and IMF), and we 
assume a constant 0.2 dex scatter in disk sizes. As we will show, our 
results are not especially sensitive to the adopted size evolution, so 
this is not a major source of uncertainty. 

2.2 "Normal" Star Formation: Relation to Galaxy Properties 

Our major assumption is the Kennicutt-Schmidt law holds at all 
redshifts, relating the (average) surface density of star formation to 
the average surface gas density 



E„ = 1.3 x lO" 4 M yr~'kpc~ 2 (- 



(3) 

,M pc-' 

with the best-fit index iik ~ 1.4 (Kennicutt 1998). The normaliza- 



tion here is corrected for our assumed IChabriei l l2003h IMF. We 
will later consider the index and normalization to be free, but for 
now take this relation as fixed. Some simple algebra shows that 
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this (assuming no dramatic evolution in disk profile shapes) can be 
written in the global form 



M yr- 



= 1.3 



(-) 
10 10 M Q J Vkpc/ 



-2(n*-l) 



(4) 



Together with the assumptions above, this defines a "steady state" 
or "quiescent" star formation rate for all disks in the model. 

We will show that the resulting SFR distributions agree well 
with those observed for normal galaxies, suggesting that these scal- 
ings are reasonable. However, we have also checked them against 
direct observations of the median SFR of disc galaxie s as a function 
of stell ar mass and redshift. The se are measured in iNoeske et all 
d2007d) from z = - 1.2 and in lPapovich et al.l d2006h at z ^2: 
comparing with the simple predictions from Equation [4] and the 
equations above yields reasonably good agreement. The combined 
dependence of / gas (Af* ) and R e (Af* ) means that Af * (Af * ) weakly 
increases with Af, in a roughly power-law like fashion, in good 
agreement with these observations; and the redshift dependence of 
/gas yields a similar increase in the normalization of the SFRs with 
redshift (and by construction, the z — normalizations are similar). 
In fact, we find that our results for the quiescent disk population 
are completely unchanged (within their uncertainties) if we simply 
adopt a parameterized fit to the observations in these papers - i.e. 
if we bypass all of the above assumptions and simply adopt a fit 
to the observed Af, (Af* | z) relations. However, this would severely 
limit the dynamic range in redshift and mass to which we could ro- 
bustly apply these models, as well as limiting the physical insight 
gained, and (most of all) would not allow for the straightforward 
predictions for merger and AGN populations. 



2.3 Merger Rates and Resulting Starburst Properties 

We next require a model for merger rates, in order to model merger- 
induced bursts of star formation and AGN activity. The methodol - 
ogy for doing so is described and tested in IHopkins et all d2009il) . 
but we briefly summarize here. 4 We assign each galaxy to a halo 
or subhalo in a simple manner f ollowing the standard halo occu- 
pation methodology described in IConrov & Wechslerl d2009l) ; en- 
suring, by construction, that the galaxy mass function and galaxy 
clustering (as a function of stellar mass, galaxy color, and physi- 
cal scale) is exactly reproduced. At a given instant, then, knowing 
the halo-halo merger rates as a function of e.g. halo mass and mass 
ratio, we can convolve this with the determined galaxy masses in 
each halo (and appropriately correct for e.g. the dynamical friction 
time delay between halo-ha lo and galaxy-galaxy merger, following 
iBovlan-Kolchin et al.ll2008l) . and obtain the galaxy-galaxy merger 
rate as a function of galaxy mass Af * , redshift z, and galaxy-galaxy 
baryonic (or stellar) mass ratio /i = Af 2 /Af 1 (defined always so that 
Af 1 > Afi, i.e. < /i < 1). The halo mass fun ctions and merger rates 
are adopted from the Milleniu m simulation 1 Springel et alj|2005bl ; 
lFakhouri&Mall2008t) : but in IHopkins et al" I d2009ih we compare 
this with a wide variety of alternative simulations and calcula- 
tions, as well as a number of differences in methodology, and show 
that these all lead to small (factor < 2) differences in the result- 
ing merger rate. In that paper, we also compare this calculation to 



4 T he approxim ate mer ger rates from the model presented 
in IHopkins et ID j2Q09ib can also be obtained as a func- 
tion of galaxy mass, merger mass ratio, and gas fraction 
from the publicly available "merger rate calculator" script at 
http : // www .cfa. harvard. edu/ ~phopkins / Site /merger calc 



a large number of observational constraints in the redshift range 
z = — 2 (see references therein), and show that the two agree well; 
adopting a parameterized fit to the observed major merger rate from 
most observations yields an identical result in our calculation here 
(but does not have the convenience of being easy to extrapolate 
to arbitrary mass ratios /i or redshifts z). Note that, again, we be- 
gin from just the star-forming galaxy luminosity function - "dry" 
mergers of quiescent systems will not produce interesting starburst 
or AGN activity. 

In a merger, gravitational torques lead to gas in the disc rapidly 
losing angular momentum to the nearby stars, and falling inwards 
( iBarnes & Hernquisd[l99ll . 1 19961) . The rapid increase in the cen- 
tral gas den sities drives a massive starburst dMihos & Hernquistl 
Il994ll996t) . Here, we assume that every merger induces a star - 
burst and corresponding AGN activity. In IHopkins et alj d2009bl) . 
as well as a number of other studies (e.g. Idi Matteo et alj 120071 ; 
ICox et al 1 l2008l ; lHopkins et alj|2009hl and references therein), the 
resulting total starburst mass/amplitude and peak QSO luminosity 
are quantified as a function of merger properties, from a suite of 
hundreds of high-resolution hydrodynamic galaxy merger simula- 
tions. These simulations span a range in the relevant properties: 
redshift, merger mass ratio, orbital parameters, galaxy structural 
properties, a nd gas fractions; and they include prescriptions based 
on the same iKennicuttl dl998l) law for dynamic star formation, as 
well a s black hole accretion and feedback from s upernovae and 
AGN dDi Matteo et~al]|2005l ; ISpringel et alj|2005ij) . Together this 
allows for a full sampling of the interesting parameter space, and a 
simple, direct parameterization of the resulting burst properties. 

Despite the complex physics involved, it is shown therein that 
the average burst scalings can be represented in analytic form, mo- 
tivated by basic gravitational physics. We adopt the full scalings de- 
rived therein, but note that the important parameter, the total mass 
of gas that loses angular momentum and participates in the cen- 
tral starburst, scales (to lowest order) with the simple relation (after 
averaging over a random distribution of orbital parameters) 



Afburst ~Af gas ft(l-/ga 8 ) 



(5) 



The scaling with merger mass ratio /i represents the declining ef- 
ficiency of angular momentum loss in more minor mergers; the 
scaling with (1 — /gas) comes from the fact that the torques that 
remove angular momentum from gas are primarily internal, from 
stars in the same galaxy - a pure gas merger would simply yield a 
new disk, not a compact starburst. 5 Adopting the simplified scaling 
above, i n fact, yields very sim ilar results to the full scaling pre- 
sented in lHopkins et alj d2009bl) for Afburst as a function of fj,, Af gas , 
/gas, and orbital parameters (we assume random orbital inclinations 
and parabolic orbits, motivated by cosmological simulations). Mo- 
tivated by the simulations (or e.g. allowing for the full distribution 



5 The physics of these scalings, particul arly that with gas fraction, is dis- 
cussed in detail in lHopkins et al] l2009bl) . In short, hydrodynamic torques 
and pressure forces are negligible, and direct torques from e.g. the sec- 
ondary galaxy and halo are suppressed by a tidal term ~ (r/R c ) 3 , the short 
time of close passage (much less than the several dynamical times needed 
to continue strong gas inflows), and the fact that they are out-of-resonance 
with the primary gas disk. Moreover these torques are just as likely to in- 
crease as decrease the gas angular momentum. As shown in that paper and 
earlier iBarnes & Hernquist 1996; Barnes|[l998l) . this means that the stellar 
disk in the same galaxy (with fractional mass (1 — /gas)), being in direct 
spatial proximity and resonance, always dominates the torques driving an- 
htgullar momentum loss in the gas. 
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of orbital parameters), we adopt a constant 0.35 dex scatter in Mburst 
at fixed galaxy properties (with of course the limit A/burst < M gas ). 

In a burst, there is some non-trivial time-dependent lightcurve 
or star formation rate versus time. Since we are considering the 
statistical distribution of luminosities, we do not need to know 
the exact time-dependent form of this function; rather, the impor- 
tant quantity is the distribution of times spent at different lumi- 



We find that it can be conveniently represented by the scaling 



nosities . Examples of this are shown in detail i n [Hopkins et al 



2006ch . and sim i lar qu antities are presented in Idi Matteo et al 



20071) : ICox et alj ( feooj) . We find that, integrated over the history 



of a burst, this function can (on average) be conveniently repre- 
sented by the simple function 



dr 



dlogM, 



: /burst In 10 exp ( - M * } , 

*- Mburst ft burst > 



(6) 



where, fitting to the simulations, we find fburst ~ 0.1 Gyr, nearly 
independent of galaxy mass and redshift. This functional form is 
characteristic of a rapid, exponential rise from low SFR to a peak 
in the burst, with a burst lifetime of order fburst- The constancy 
and normalization of this lifetime is a simple consequence of the 
observed dynamical times in the central regions of galaxies, and 
the fact that these dynamic al times scale weakly or not at all with 
mass and redshift (se e e.g. iBell & de Jongll200ll : lMcGaugh|[2005l ; 
ICourteau et alj|2007l) . Given some merger rate, and correspond- 
ing rate of "creation" of bursts of a given mass (dn(Mburst)/df), 
the observed number density of bursts at a given M„ is simply 
given by the convolution of this rate with the lifetime above, i.e. 
An (Mburst )/d log M* = (dw(Mbu fs t) /At ) (df /dlogM*) (for m ore ex- 
plicit details of this methodology, see H opkins et alj2006ah . 

Note that th e numbers above are somewhat different from 
those presented in lHopkins et al. I J2006d) . However, in that paper, 
we were considering the total distribution of star formation rates 
that would be observed over the entire duration in which a system 
might be identified as a merger or interacting pair. The lifetime of 
that phase is much longer, ~ 1 — 2 Gyr, and (by time and by to- 
tal mass) most of the star formation comes from the "normal" star 
formation that would be associated with the two merging dis ks in- 
dependently (see e.g. Idi Matteo et~ai]|2007l ; ICox et alj|2008h . De- 
pending on the observational criteria used to identify mergers, of 
course, this definition of star formation "in mergers" may be of in- 
terest. It is, however, a subset of the "quiescent" star formation for 
the most part, and is distinct physically (and very distinct in terms 
of the imprint that it leaves on galaxy stellar populations, kinemat- 
ics, and structural properties) from the short lived, compact burst 
specifically induced by the merger. 



2.4 AGN and Quasars 

Given some merger, quasar activity is also excited; to lowest order 
in simulations with AGN-feedback, the peak bolometric luminos- 
ity of the AGN is tightly coupled to the total bulge mass that will 
be formed from a disk-disk merger. The total bulge mass is the 
burst mass (discussed above), plus the violently relaxed stellar disk 
mass, which is simply M re iaxed ~ fiM, in simulations and from sim- 
ple gra vitational physics considerations (again, see IHopkins et al .1 
l2009bl) . At fixed bulge mass, the peak AGN luminosity, corre- 
sponding to the Eddington limit of the maximum BH mass, is cou- 
pled to this bulge mass as it mus t overcome its binding e nergy in or- 
der to halt continued growth. In lHopkins et al. U2007bfi) . we show 
how this scales in simulations with bulge mass and other properties. 



Lpeak,QSO~4.6X 1O"L0 X (1+z) 



0.5 / Mburst + Mrelaxcd N 



10 10 M, 



© 



(7) 



The latter scaling simply reflects the fact that the peak/final BH 
mass scales roughly linearly with total bulge mass, in both obser- 
vations and simulati ons. The (1 +z)°' 5 scali ng comes from the sim- 
ulations discussed in lHopkinsetal1 ( l2007ah : it comes from the fact 
that galaxies at high redshift, being both more gas-rich and more 
compact, require more "work" to be done by the AGN before it 
can self-regulate its luminosity/BH mass, and so yield higher BH 
masses at otherwise fixed bulge mass. We refer to that paper for 
more details, but note that the other parameterizations of this evo- 
lution (discussed therein) yield nearly identical results. Likewise, 
other models for AGN self-regulation at high masses and/or lumi- 
nosities predict a similar maximum BH mass as a function of host 
galax y properties dSilk & Reesj|l998l : lMurrav et alj|2005l ; IShankaj 
|2009|) . and the resulting "cutoff" in the AGN luminosities at high 
masses (owing to self-regulation combined with a cutoff in the 
depth of host galaxy poten tial wells) is similarly predicted in e.g. 
iNataraian & Treisterl ll2009). Moreover, such moderate evolution in 
Mbh/M« is suggested by a number of observations (e.g.lPeng et al] 
Dd lWoo et al.l2 



120061 : IWoo et alhoOol ; iTreu et alj2007l ; ISalviander et al.l200l and 

references therein). Motivated by the simulations and the observed 
BH-host correlations, we assume a constant 0.3 dex scatter in these 
relationships. 

In what follows, we consider only AGN induced in merg- 
ers, but stress that this does include non-trivial contributions from 
minor mergers down to e.g. mass ratios of ~1:10. There is con- 
siderable debate regarding whether or not entirely non-merger 
processes such as stellar b ars (e.g. IShlosman et al 1 ll989l ; |jogeel 
120061 : 1 Younger et al.l2008al and refere nces therein) and/or stochas- 
tic encounters with molecular clouds faopkins & Hemquistll2006l: 
iNavakshin & King! |2007|) might drive significant AGN activity. 
However, it is generally clear both from observations and from 
simple theoretical considerations th at the resulting AGN would be 
important only at low luminosities. iHopkins & Hernquistl j2009al) 
compile both empirical and theoretical estimates of the lumi- 
nosities below which non-merger processes dominate AGN fu- 
eling and find consistently that this is excl usively in the tradi- 



tional Seyfert regime (L bo i < 10 12 L^; see llvlalkan et al] 1 19981: 



ICanalizo & Stockton! 1200 ll; iDunlop et al] 120031; iKauffmann et al 
| 2003| : iFlovd et al] |2004 iHutchings et al] 120061: IZakamska et al 
20061 (20081 ; iRigbv et all (20061 : iGuvon et al] I2OO6I: lUrrutia et al 
[2008J). This is clear from integral constraints; bulges formed in 
bars ("pseudobulges") or non-merging bulges with bars or central 
molecular gas concentrations domi nate only at low galaxy masses , 
in galaxy types of Sb/c and later |Kor mendv & Kennicutt |2004|; 



Fishdl2006l : lAllen et al.|[2006l ; iDriver et al] 2007; Fis her & Drorvl 
2008b . Given the observed BH-host correlations, this corresponds 



to BHs with masses < 1O 7 M , or maximum luminosities at the 
Eddington limit of Lboi = 3x 10" Lq. But as we will show, AGN 
are significant in the IR luminosity function only at the highest 
luminosities, L > 10 13 Lq - i.e. BHs with ~ 10 9 M Q at Edding- 
ton, with accretion rates of > IOMq yr~' . Since it is unlikely that 
these extreme systems are powered in non-violent events, our ne- 
glect of non-merger induced AGN makes little difference to our 
predictions. We have, in fact, explicitly checked whether includ- 
ing them (according to the m odel luminosity functions predicted in 
IHopkins & Hernquistll2006T) makes any difference, and find it only 
increases the very low-luminosity contributions of AGN by a factor 
of ~ 2 — 3, far less than the ~ 3 — 4 orders of magnitude required 
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to substantially change our conclusions. Likewise, simply adopt- 
ing the observed AGN bol ometric luminosity fun ctio ns - including 
all obs erved AGN - from IHopkins et al ] d2007ch or Shankar et al. 

with our estimated template spectra and obscured fractions, 
we find no significant difference. 

For a given peak BH mass or peak luminosity, we simply re- 
quire again the distribution of luminosities corresponding to the 
average lightcurve, in order to construct the number density as a 
function of luminosity. These lightcurves and the resulting distri- 
bution of time spent at different AGN luminosities (both bolomet- 
ric and in various ob served bands) have be en e xtensiv ely discussed 
in a series of papers dHopkins et al . 2005a d U ljjOol lbf). We adopt 
the Schechter-function parameterization therein, 

* - 0.22 Gyr ( I**)" exp { - (***) \ 

d log Lbol V Lpeak / *■ V Lpeak / > 

a « -0.44 + 0.21 log(L pMk /lO I2 L ) (8) 

Again, this is taken from simulations, but in those papers, it is 
shown that this yields very good agreement with the observed dis- 
tribution of AGN Eddington ratios, host masses and luminosities, 
and the evolution of the AGN luminosity function (most recently , 
seelVolonteri et alj|200d iForeman etal . 2008; Bonoli et al.l l2008l; 
Marulli et alJ 120091 : iHopkins et ail l2008dl : iHopkins & Hemquistl 
2009al lbl). 

2.5 Construction of IR Luminosity Functions 

Finally, given these predicted SFR and AGN bolometric luminosity 
distributions, we need to convert to the observable quantity, namely 
total infrared luminosity. Because we are not attempting to model 
the full SEDs and dust physics of these systems, the only quantity 
that we can robustly predict is the total infrared luminosity, Lir, de- 
fined as the integrated luminosity from 8 — 1000 /jm. In the case of 
SFR d istributions, we adopt the s imple conversio n from lKennicuttl 
dl998h , corrected for our adopted [chabrierl d2003l) IMF, of 

L^l.lxlO^^). (9) 

Note that more sophisticated (e.g. luminosity-dependent) conver- 
sions have been proposed, but since this choice is used to calibrate 
the gas surface density-SFR surface density relation, we adopt it 
for consistency. In any case, alternative formulations largely devi- 
ate from the above only in non-starburst or lower IR-luminosity 
galaxies, where absorption is weaker, but these are not particularly 
i mportant for our compariso ns here, and experiment ing with those 
in lBuat & Burgarellal Jl998h and ljonsson et al.l J2006l) yields almost 
no difference at Lir > 10" Lq. 

An advantage of our semi-empirical model is that our conclu- 
sions do not depend significantly on the adopte d stellar IMF, fo r 
typical c hoices . Altering the IMF bet ween e.g. IChabrierl d2003l) . 
iKroupal d2002h . IScald dl986h or e.g. ISalpeterl dl955h generally 
amounts to systematic changes in the mass-to-light ratio M«/L by 
up to 0.3 dex. However, because we begin with observed galaxy 
properties and calculate observed luminosity functions, this sys- 
temati c dependence can cels out. Specifically, in adopting the ob- 
served [Kennjcutjl d 19981) relation, the observable quantity is the lu- 
minosity surface density - this factor enters in the conversion to a 
SFR surface density. But then, converting the resulting SFR to an 
observed luminosity, the same factor enters, cancelling out. Explic- 
itly re-calculating our predictions with other IMF choices confirms 
this. The only residual effects are second order, and relatively weak 



- for example, changing the implied gas exhaustion timescale leads 
to slightly different dynamics of the gas on small scales in mergers. 
Such details are outside the scope of our comparison here, and in 
any case amount to smaller effects than our systematic uncertain- 
ties. In an a priori model for star formation, on the other hand, the 
factor would enter fully. The IMF will only present a systematic 
source of uncertainty in our predictions if it evolves significantly 
with redshift or galaxy properties, or if it is extremely top-heavy, 
possibilities we discuss further below. 

For AGN, the conversion from bolometric to IR luminosi- 
ties is somewhat more complex - unobscured (Type 1) AGN re- 
radiate only a fraction ~ 1/40— 1/20 of their bolometric lumi- 
nosity in the FIR, and are thus negligible for the luminosity func- 
tions here. We adopt the empirically ca lculated obscured fraction 
as a function of quasar luminosity from lGilli et alJ J2007t) . and as- 
sume that the obscured bolometric luminosity is re-radiated in the 
IR; this allows us to convert our predicted bolometric QLF to an 
IR QLF of obscured quasars. Technically, not all of the luminos- 
ity will be obscured, of course, but we find that e.g. using the 
full distribu tion of column den sities as a function of quasar lumi- 
nosity from lUeda et alJ d2003l) to atte nuate a template AGN S EP 
yields a very similar answer (see also lFranceschini et "al]|2005h . as 
does using a mean X-ray to IR bolometric correction of obscure d 
AGN dElvis et al.ll994l ; IZakamska et al.l2004l ; |Polletta et alj2006h . 
The obscured A GN fraction at high luminosities remains uncertain; 
lHasingerl d2008l) argue tha t it could be lower by a factor of several 
than the lGilli et alJ J2007h estimate at the highest luminosities (al- 
though this is redshi ft-dependent and at z > 1.5, the two estimates 
agree well), whereas [Daddi et all d2007l) argue that the number of 
obscured quasars should be a factor of ~ 2 — 3 higher. These un- 
certainties are generally comparable to (or up to factor ~ 2 larger 
than) the uncertainties from the choice of stellar mass function, and 
we discuss some implications below. 



3 STAR FORMATION RATE DISTRIBUTIONS AND 
INFRARED LUMINOSITY FUNCTIONS 

3.1 Basic Predictions 

Figure[2]shows the resulting predicted total-IR luminosity functions 
from z = — 6, divided into the contributions from "normal" (non- 
merging, quasi-steady-state) star-forming systems, merger-induced 
"bursts," and obscured AGN. As discussed in § [2] we have re- 
calculated our model adopting at least two different stellar MF de- 
terminations at each redshift; the range between the two at each 
redshift is shown by the shaded range, and is representative of the 
scatter in different observational estimates. Unsurprisingly, this un- 
certainty is substantial at high redshifts. We also add (in quadra- 
ture) a systematic factor ~ 2 uncertainty in galaxy-galaxy merger 
rates, representative of the systematic theoretical an d observational 
uncerta inties as estimated from the compilations in IHopkins et al.l 
d2009il) . We add a factor 1 .5 additional uncertainty in the AGN ob- 
scured fractions, aga in representative of systemat i c observational 
uncertainties (see e.g. IHopkins et alj2005dl2009dl ; lshi et al.l20o3 ; 
iTreister & Urry||2"006l ; lHasingerll2008l , and references therein). We 
also note that direct observational constraints used for our mod- 
els of the galaxy stellar mass function are either non-existent or 
extremely uncertain above z > 4; we extrapolate the fitted LF pa- 
rameters from z = 2 — 4 into this redshift range, and so the resulting 
predictions should be treated with the appropriate caution. 

We compare with observations of the IR luminosity functions 
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Figure 2. Total (8 — 1000 /^m) IR luminosity functions as a function of redshift. We show the model contribution from "normal" (non-merging) star-forming 
disks (green), merger-induced starbursts (blue), and obscured AGN (red). The range in the total (summed) LF is shown with dotted black lines. Shaded 
ranges reflect the uncertainty from different stellar mass function observations used in constructing the model. Points show observational estimates from 



Saunders et al.l (19901 magenta stars), ISoifer & NeugebaueJ jl99lL blue triangles), [Yun^ni^ bOOll blac k circles), iHuang et all tool magenta squares), 
violet diamonds) jQrpu^TetaiT i l2007l orange inverted triangles uMagnelli et alj|2009l black *'s), BabbedgTetalTl200d . red pentagons), 



Le Floc'h et al 



Chapman et al 



20051 dark green +'s). 



where available, from z = — 3. Note that all of these are corrected 
to a total IR luminosity fro m observations in some band; we adopt 
the corrections compiled in IValiante et all J200gj> . but emphasize 
that some caution, and at least a systematic factor ~ 2 uncertainty in 
Lir, should be considered in estimates from most if not all observed 
wavelengths. The agreement between the total predicted LF and 
the observations is generally reasonable, at most redshifts. At the 
highest luminosities a nd redshifts, specifical ly the sub-millimeter 
population observed in lChapman et alj J2005t> . we appear to under- 
predict the abu ndance of bright systems , but these observations are 
very uncertain. I Austermann et al. (2009]), for example, find that the 
millimeter number counts in these surveys are strongly affected by 
cosmic variance, and may be factors of several larger than the cos- 
mic mean. We discuss this in § [4] but note for now that these sys- 
tems contribute relatively little to the global SFR density at these 
redshifts. 

At all redshifts, "quiescent" galaxies dominate the LF at low 
luminosities and high space densities, reflecting the abundance of 
star-forming disks and relative rarity of mergers and quasars. At 
higher luminosities, eventually merger-induced star formation and 
AGN activity become dominant, as expected in order to explain the 
most extreme (but short-lived) bursts of star formation. However, 



both types of systems increase in luminosity with redshift from z = 
— 3 in similar fashion. 

As discussed in § 12.51 the obscured AGN fraction is some- 
what uncertain. However, Figure [2] shows that even at extreme lu- 
minosities, the contribution from obscured AGN is comparable to 
that from merger-induced starbursts. Thus, in terms of the total IR 
luminosity function, even an obscured AGN fraction of zero would 
only lead to factor ~ 2 changes in the predicted bright-end num- 
ber densities (smaller than the uncertainties owing to the choice of 
mass function, for example). 

To facilitate future comparisons with observations, we present 
the corresponding predictions in Appendix lAlfor the IR luminosity 
function in various specific rest-frame wavelengths. However, we 
stress that these are not direct predictions of the model - a proper 
model for the SEDs will depend on full radiative transfer models, 
applied to the simulations as a function of time and galaxy prop- 
erties (these will be presented in future work). Here, we simply 
convert total IR luminosities to wavelength-dependent luminosities 
using the same bolometric corrections used to convert the observa- 
tions in Figure[2] 
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Figure 3. Best-fit parameters for the luminosity functions in Figure[2] given 
a double-power law formulation (Equation HQ) . Points show the best fit 
(with uncertainty reflecting the allowed range in Figure |2) at several red- 
shifts, lines the overall maximum likelihood fits from Table[T] Parameter de- 
generacies are such that the best-fit curves can be systematically slightly off- 
set from the fits at each redshift. Top: Break luminosity. Rising gas fractions 
drive the increase; but for mergers, the results asymptote to a maximum ow- 
ing to the physics of very gas-rich mergers (see text). Second from Top: Nor- 
malization. Modulo a decrease at the lowest redshifts, this is approximately 
redshift-independent. Second from Bottom: Faint-end slope. Again, the red- 
shift dependence is weak. The apparent large evolution for star-forming sys- 
tems is somewhat degenerate with the evolution in the bright-end slope - a 
fit where both are held constant is, in fact, acceptable. Bottom: Bright-end 
slope. Again, relatively flat with redshift (per note above). Star-forming sys- 
tems fall off in number density at high luminosities much more steeply than 
mergers or quasars, reflecting the exponential cutoff in the mass function. 



4>: 



dn 



dlogL {L/L*)<* + {L/L*Y 



(10) 



where the parameters c/>, (normalization), L* (break luminosity), a 
(faint-end slope, i.e. $ oc L~ a for i<L»), and j3 (bright-end slope, 
i.e. $ oc L~" for L>L,) depend on redshift, with that dependence 
conveniently approximated as 



logL* 


= L + L'i + L"i 2 


log</>* 




a 


= ao + a ? + a £ 


P 


= /3o + f3'Z + f3"f 


€ 


slog(l+z) 



(11) 

(Note that log here and throughout refers to log 10 .) We perform this 
fit using only our results up to redshift z = 4, as the HOD con- 
straints used to build the model have to be extrapolated at higher 
redshifts. In TableQ] we quantify the uncertainty in each parameter; 
this reflects the systematic theoretical uncertainties shown in Fig- 
ure [2] (the shaded range), with the appropriate covariance between 
parameters taken into account (for this reason, fitting the redshift 
evolution with free parameters up to second-order in £ leads to rel- 
atively large uncertainties in the fit results). We also illustrate the 
best-fit parameters as a function of redshift in Figure[3] 

The behavior seen in each parameter reflects that discussed 
above; the bright and faint-end slopes, and normalization c/>„, 
evolve relatively weakly with redshift. 6 In fact, we can find rea- 
sonable fits within the theoretical uncertainties that hold these pa- 
rameters fixed with redshift. But the break luminosity L, evolves 
rapidly, as cc (1 +z) for z < 2 in all populations, then levels out 
to a maximum at higher redshifts. We do see this flattening from 
z ~ 2 — 4, hence the quadratic term here; although given our z < 4 
limit (log ( 1 + z) < 0.7), we are only just sensitive to the quadratic 
terms in £ (and see no significance fitting higher-order terms). 

At all redshifts, L» is higher for merger/AGN populations 
(> 10 12 Lq) relative to normal galaxies; but the space density c/>„ 
is much lower (by a factor of ~ 100 — 300). The bright-end slope 
of the normal population is steep, reflecting the rapid exponential 
cutoff in the galaxy mass functions; the bright-end slope in the 
merger/AGN populations is much more shallow, ~ 2.5 — 3 - such 
a slope is, in fact, very simila r to the observed bright-end slope of 
the b r ightest IR population s dSanders et all 1988a: Sa unders et al.l 
ll990l : IChapman et alj|2005l) and to the well-constrained bright-end 
slope of the quasar lum i nosity functions from redshifts z ~ — 6 



Fan et al. 2004; Brown et al.l 2006; Richards et al]|2006bl: 



iHopkins et al]|2007d ; IShankar et ai]|2009l : ICroom et ai]|2009l and 
references therein). 

For comparison, we also consider fitting the LFs to a modified 
Schechter function parameterization, namely 



<1> 



(12) 



with the same assumed form for the evolution in the fit parameters 
with redshift. This is akin to a standard Schechter function except 



3.2 Fitting Functions to the Predicted LFs 

For the sake of comparison with future observations, we provide 
fits to the model predictions in Table Q] We find that the predicted 
IR LFs for each type can be reasonably represented by a double 



6 The apparent "jump" in 0* at z ~ 0.3 owes partly to real evolution in the 
observed input mass functions, but mostly to parameter covariance (here 
between tf>„, L», and a). Accounting for this covariance, the change in c/>* 
from z = 0.2 — 0.3 is only significant at 1.5 — 2cr, and a smoothly evolving 
<f>* provides just as good a fit. 
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with the addition of a bright end "slope" term /3, where /3 < 1 al- 
lows for a less-steep falloff at high-L than would be predicted by 
a standard Schechter function (/3 = 1). We find that, because the 
functions shown in Figure|2]do not have sharp "breaks" character- 
istic of a double power-law, this provides a marginally more ac- 
curate representation of the LF shape. However, the difference is 
small, and direct interpretation of the parameters in Equation 1121 
is complicated by serious fitting degeneracies. With this choice of 
functional form, we find the second-order redshift evolution terms 
make little difference to the fits, and so - given the steep parameter 
degeneracies involved - do not free the higher-order terms in the fit. 
Interestingly, the total luminosity function obtained by summing 
the contributions from each component is better represented with a 
double power-law, as opposed to the modified Schechter function. 



3.3 The Luminosities & Space Densities of Population 
Transitions 

We explicitly quantify the "transition point" between the domi- 
nance of one population or another as a function of redshift in 
Figure [4] Specifically, we define this as the point where the lumi- 
nosity functions from different populations in Figure [2] cross. For 
example, the transition luminosity or space density between dom- 
inance by normal disks and mergers is given by the point in Fig- 
ure [2] where (Lir | normal) = 0(Lm | burst); above this Lir (and 
below the corresponding 0(Lir)), 0(Lir| burst) > tf> (Lir | normal), 
at lower luminosities and higher space densities the opposite is true. 
Likewise, we can define the transition luminosity or space density 
where obscured AGN become numerous than star-formation domi- 
nated systems, <?!>(Lir | AGN) = ^>(Lir | normal) + <f> (Lir j burst). The 
uncertainties in Figure [2] are translated to corresponding uncertain- 
ties here. 

Our comparisons generally affirm the conventional wisdom: 
at low redshift, mergers dominate the ULIRG and much of the 
LIRG populations, above a luminosity ~ 10 1L5 Lq. Heavily ob- 
scured (potentially Compton-thick) AGN (in starburst nuclei) be- 
come a substantial contributor to IR luminous populations in the 
most extreme > a few x 10 i2 Lq systems (nearing hyper-LIRG 

> 10 13 Lq luminosities which are common bolometric luminosities 
for > 10 Mq BHs near Eddington, but would imply potentially 
unphysical > 1000 Mgyr~' SFRs). At higher redshifts, disks are 
more gas-rich, and thus have characteristically larger star formation 
rates, dominating the IR LFs at higher luminosities. By z ~ 1, most 
LIRGs are quiescent systems, and by z ~ 2, only extreme systems 

> a few x 10 12 Lq are predominantly mergers/ AGN. 

This appears to agree well with recent observations. First, 
consider the results of systematic morphological studies of IR- 
bright sources as a functi on of their luminosities, at low redshifts 
dSanders & Mirabel 199^ . which affirms the conclusion that - lo- 
cally - the brightest LIRGs and essentially all ULIRGs are merging 
systems, while less-luminous systems are not (see also references 
in § [TJ. At high redshifts, sim ilar studies have now been performed 
(see e.g. lTacconi et alj[2008l , and references therein). They too find 
that the brightest sources are almost exclusively mergers, but with a 
transition point (from non-merger to merger-dominated) an order- 
of-magnitude larger in luminosity. Other morphological studies at 
inter mediate redshifts z ~ 0.4 —1.4 have reached similar conclu- 
sions ( iBridge et ai]|2007l) . 

Other studies have attempted to separate the contributions 
of star formation and (obscured) AGN. At low redshifts, we 
find similar r esults from observational comparison of emi ssion 
line strengths dSanders & Mirabej [19961 ; iKewlev et al.ll2006l) . ob- 



servations of the strength of observed PAH features I Lutz et al 



1998| ; IVeilleux et alj2009bh . full SED template fitting dFarrah et al 



20031), or indirect comparis on with Type 2 AGN luminosity func- 
tions dCharv & ElbazkOOll) . In each case, these studies find that lo- 
cal "normal" ULIRGs are star-formation dominated, but extremely 
rare systems approaching Hyper-LIRG luminosities tend to be 
AGN-dominated. At high redshifts, there have recently been at- 
tempts to apply similar methodologies, especially co mparison of 
PAH s trengths, and we show such an estimate from ISaiina et al.l 
d2007t) . who find a similar transition from star formation to AGN 
do minance as at low r edshift, at a factor of several higher luminos- 
ity. lYang et al] d2007t) measure dust temperature distributions and 
positions on the far IR-radio correlation for a sample of ULIRGs 
over the redshift range z = 0.3 — 1; they find that below Lir = 
10 12 ' 4 Lq (Lfir = 10 12 ' 25 Lq), the systems appear star-formation 
dominated, while above Lir = 10 12 9 Lq, the IR luminosities are 
dominated by A GN. At z = 2, the same constrain t s suppo rt the 
conclusion from lSaiina et alj d2007l) ; lYounger et alj d2009K ) show 
that samples of ~ 2 - 8 x 10 12 Lq ULIRGs at z = 2 follow the 
local far I R-radio correlati o n, indi cating they are starburst domi- 
nated, but iBussmann et alj d2009al lbl) find that by luminosities of 
~2x 10 13 Lq, IR samples are dominated by warm dust sources 
more likely to be (post-merger) AGN. Of course, changing the as- 
sumed number of obscured AGN, as a function of luminosity or 
redshift, will correspondingly shift the predicted transition point; 
the agreement seen here suggests that the correct number is prob- 
ably not very different from that adopted here. For fur t her details, 
we ref er to the above as well as IChapman etal.ld2005l) . lDevetal.l 
d2008l) , and lCasevetai]d2009l) . 

The transition point between non-merger and merger domi- 
nance of the luminosity function shifts to larger luminosities at high 
redshifts, even though gas-rich merger rates increase rapidly. The 
evolution in gas fractions, which drives up both disk star formation 
rates and merger-induced bursts similarly, is the dominant effect; 
the evolution in merger rates is also not .so rapid as to dominate 
the population at redshifts z < 2. Moreover, as noted in § [2] al- 
though disk star formation rates (at otherwise fixed properties) in- 
crease monotonically with their gas fractions and hence gas surface 
densities, merger-induced bursts and quasar episodes can decline 
in efficiency in extremely gas-rich systems, because the gravita- 
tional torques that allow for such bur sts depend on a sizab le dis- 
sipationless (stellar) disk component dHopkins et al.l|20 09b). As a 
result, merger-induced bursts do not grow in importance as rapidly 
as might naively be expected from analysis of e.g. halo-halo merger 
rates. 



3.4 Corresponding SFRs and Bolometric Luminosities 

In Figure [5] we reproduce our model from Figure [2] but show in- 
stead the distribution of star formation rates. Of course, AGN are 
not present here, since although there will be star formation in their 
hosts, the AGN IR luminosity itself is not from star formation. 

Figure|6]shows the corresponding bolometric luminosity func- 
tions. For the star-forming systems, this is essentially identical to 
the total-IR luminosity functions, as the total IR emission domi- 
nates the bolometric luminosity in at least the luminous, IR-bright 
end of the distribution (of interest here). However, for the AGN, the 
bolometric emission is considerably larger than the IR emission. 
We therefore highlight the AGN predictions. We compare these 
to the large compilation of obser vations used to derive bolomet- 
ric quasar luminosity functions in lHopkins et al. | d2007cb (see ref- 
erences therein). Similar results have been obtained in other compi- 
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Figure 4. Top Left: Total IR luminosity threshold above which the predicted IR luminosity functions in Figure [2] transition from being dominated by non- 
m erging ("quiescent") disks to merger-i nduced star formatio n/bursts. We compare observational estimates from morphological studies of bright-IR sources 



t") 

in Sanders & Mir abel H99(l circle) and lTacconi et alj fc008t triangle). Bottom Left: Same, but showing the space density threshold (<!>, in Mpc -3 log -1 Lir) 
of the same transition. Top Right: Luminosity threshold above which the IR luminosity functions transition from being dominated by star for mation to being 
domin ated by AGN. Points show the ob served estimates from comparison of PAH feature strengths and emission line strength template fittin g in | C harv & Elbaj 
1200 ll square). Fveilleux et al l 1 2009b a, diamond), and lSaiina et alj 120071 star), and from comparison of the far IR-radio correlation in lYang et*alTl2007h . 
Bottom Right: Same, in terms of the space density threshold. As gas fractions increase with redshift, star formation rates increase in all systems. As a result, 
the threshold for merger-dominance grows from bright LIRGs at z ~ 0, to bright ULIRGs at z ~ 1 — 2, to HyLIRGs at z > 2. In terms of space density, this 
transition is relatively constant over this range at 10~ 6 — 10 — 5 Mpc -3 log - 1 Lir (at higher redshifts, the space densities of all massive systems drop rapidly). 
At all redshifts, bright HyLIRGs (Lir > 10 13 Lq ; $ ~ 10~ 7 Mpc -3 log - 1 Lir at z ~ - 4) have a non-negligible AGN contribution. 
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Figure 5. As Figure [2] but showing the distribution of star formation rates 
at each redshift (hence no AGN contribution). 



lations l lShankar et aT]|2009l) , or from hard X-ray lumin osity func- 
tions with appropriate bo lometric corrections (see e.g. lAird et all 
2()0'J:lYLMicho et al.ll2009h . The observationally estimated bolomet- 
ric QLF agrees reasonably well with our predictions. The under- 
prediction of very low-luminosity AGN owes to our neglect of 
non-merger induced AGN (§ I2.4t ; it is clear here that these have 
a negligible impact on our conclusions. The model may also some- 
what under-predict the number density of the most luminous sys- 
tems (L boi > 10 14 Lq); this is discussed in detail in lHopkins et alj 
d2008dl) . but is sensitive to the assumed scatter in bolometric cor- 
rections, and to the existence of even a small lensed or beamed 
QSO population. 

3.5 The Effects of Different Model Assumptions 

We briefly outline the effects of several important components in 
the models adopted. For further details, see Appendix \E\ where 
we reproduce our model from Figure [2] explicitly, with different 
changes (discussed here) to the model. 
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Figure 6. As Figure [2] but in terms of bolometric luminosity (similar to 
the IR for star-forming systems, but significantly larger for AGN, which 
are highlighted here). Black points show the compilation o f observational 
data us ed to derive bolometric AGN luminosity functions in lHopkins et alj 
these should be compared to the predicted AGN LF. 



If we do not allow galaxy gas fractions to evolve with redshift 
(i.e. adopt the z — value at all redshifts), this leads to a substantial 
under-prediction of the luminosities of "quiescent" galaxies at z > 
1. In short, the existence of apparently "normal" galaxies at high 
redshifts, with ULIRG-level luminosities, requi res very high gas 
surface densities (relative to those at z = 0) if the lKennicutj il998h 
relation is to hold in some form. 

If we do not allow for disk sizes to be more compact at high 
redshift, this yields lower surface densities, and hence somewhat 
lower SFRs, but the effect is relatively minor. Because the ob- 
served size evolution (of star-forming galaxies) is weak (R,, (M„ ) oc 
(1 +z)~' 0-0 6 ') and the size evolution (at otherwise fixed proper- 
ties) only enters into the SFR at sub-linear order (for E* oc Eg a t, 
this yields M„ oc i? -0,8 at otherwise fixed properties), the total dif- 
ference is relatively small (factor ~ 2) in luminosity, comparable to 
many of the other uncertainties involved. 

If we do not allow for scatter in any quantities (e.g. disk 
sizes, gas fractions, burst masses, quasar bolometric corrections, 
and SFRs at otherwise fixed properties) - i.e. force all values to ex- 
actly trace the medians given in §|2]- this has the expected effect, 
that the rare, high-L population is significantly suppressed. These 
objects depend on the existence of some systems with relatively 
high gas fractions at high masses and high Lir relative to many of 
their other properties. 

We can also re-construct our model pred i ctions , but adopt a 
more steep power-law index for the lKennicuttl i ll 9981) relation, for 
example E, suggested by som e recent observations of 

high-redshift systems dBouche et alj|2007r ). In order to avoid over- 
producing local star formation rates (and indeed the LFs at all 
luminosities and redshifts), it is necessary to correspondingly re- 
normalize the relation: we do so such that a Milky- Way like disk, 
with effective gas surface density ~ 3 x 10 8 MQkpc~ 2 (10% gas 
fr action) has the sa me SFR as that expected from the relation fit 
bv lKennicuttl jl998T) . This amounts to a factor ~ 3. 1 lower normal- 
ization in Equation[3] with the steeper hk. Considering simulations 
with such a steeper index, the resulting burst properties are qualita- 
tively similar, but the burst timescale in Equation|6]is shorter by a 
factor ~ 2, fbui-st ~ 0.4 x 10 s yr. The AGN properties are relatively 
unchanged. Together, the results from this revised model are simi- 
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Figure 7. Distribution of SFR in galaxies of fixed mass (M* = 10" Mq) 
at z = and z = 2. We show the total (black), contribution from "normal" 
systems (green), and contribution from merger-induced bursts (blue). The 
"normal" systems here have scatter reflecting that in their radii and gas frac- 
tions (assumed lognormal). Mergers, even at z = 2, do not dominate the 
population or scatter in SFR(M»), except in the extreme wings, consistent 
with observations indicati ng small deviations about the median M„ -SFR 
trend in normal galaxies (Noeske et al. 2007a). Observations probing be- 
tween ~ 2 — 3 a in the wings are needed to see the merger "tail." 

lar to our default model - however, the steeper index leads to more 
star formation in the very high gas density systems at high redshift. 
This actually somewhat improves the agreement with the observed 
number densities of the most luminous systems; however, the dif- 
ference is ultimately within the range of our other uncertainties, in 
particular the number density of the most massive galaxies. 

Again, to facilitate future comparisons, we provide in Table [T] 
fits to the same double-power law functional form for the predicted 
luminosity functions in both the case of no gas fraction evolution, 
and the case of a steeper Kennicutt-Schmidt index. 

3.6 The Luminosity/SFR Density: Contribution of Mergers 
and AGN 

3.6.1 Luminosity Densities: Predictions 

In Figure [7] we examine the distribution of SFRs at fixed galaxy 
stellar mass (for an ~ L», M* = 10" Mq system), for the standard 
model used in constructing Figure H at z = and z = 2. We sep- 
arately show the distribution of SFRs from the "quiescent" (non- 
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merger) systems at that mass, and from the merger-induced bursts. 
The scatter in non-merger systems comes from the distribution of 
gas fractions and effective radii, at a given stellar mass and redshift. 
Obviously, SFRs are systematically higher at z = 2, and the merger 
contribution is relatively larger, as merger fractions observed have 
increased from ~ 1% at z = to ~ 10% at z = 2 (see e.g . 
Bundv et al J 120091: IConselice et alj 120091; iKartaltepe et al J 12001 
Lin et alj2008l;lBluck et alj2008l ; lHopkins et al.l2009il ; |jogee etal] 



2009l : lBridge et al.||2009l. and references therein). However, at both 
redshifts, the merger contribution is relatively small, and although 
it dominates the tail at very high SFR at fixed mass, it constitutes 
much less than the > 30% of the population needed for it to bias the 
~ lcr scatter in SFR(M»). Various observations have shown that 
there is a tigh t sequence of SFR w ith galaxy mass in star-forming 
systems (e.g. iNoeske et al. I l2007al) . with small scatter < 0.3 dex, 
similar to that predicted here. This presents a constraint on the role 
of merger-induced bursts in affecting SFRs, but one easily satisfied 
here - far from affecting the scatter at the 1 a level, one has to ob- 
serve the scatter at a level between 2 — 3 a in the high-SFR "wings" 
of the distribution at fixed stellar mass before the merger-induced 
tail would be evident. 

Figure [8] combines the LFs predicted in Figure [2] to show the 
total infrared luminosity density, and corresponding total SFR den- 
sity, of the Universe as a function of redshift. Approximate fits to 
these predictions can be obtained by simply integrating the fitted 
LFs in Table [T] The luminosity density in quiescent systems dom- 
inates the global total at all redshifts. Merger-induced bursts con- 
tribute a relatively small fraction to the global SFR density; ris- 
ing from ~ 1 — 5% at low redshifts z ~ to a roughly constant 
~ 4- 10% atz > 1. 

The contribution from obscured AGN is at most comparable to 
that from merger-induced bursts, and in general a factor of ~ 2 — 3 
lower (of course, the conversion to SFR density is not valid for 
AGN, as the IR emission is powered by accretion; they should be 
compared to the total IR luminosity density only). It is unlikely 
that obscured AGN contribute more than ~ 5% to the global IR 
luminosity density, even assuming a generous near-isotropically 
obscured fraction of ~ 1/2 at high luminosities (large given the 
observational constraints from e.g. iGilli et al .1 12007k iTaier et al.1 
2007l:lDaddi et alJ^OOTllHickox et alj2007l:ICaccianiga et alj2008l; 
Treister et al. 2008; Menendez-Delmestre et alj2009l;lTreister et all 



20091 ; Malizia et al. 2009; Trichas et al. 2009). In fact, allow 



ing the entire bolometric AGN luminos i ty den sity estimated in 
IHopkins etal] j2007ch and lShankar et all J2009h to be re-radiated 
in the IR increases the contribution of AGN by only a factor ~ 3. 
Even under conservative assumptions, then, the contribution from 
obscured AGN is much less than the other current statistical and 
systematic errors in the estimation of the global SFR density, and 
infrared-derived SFR densities are not likely to be significantly con- 
taminated by AGN. This question has been studied via other means, 
as well - for example, in X-ra y background synthesis models - with 
simil ar conclusions (see e.g. iTreister & Urry||2006l : iTreister et al.1 
l2008h . 

These are global statements - at a given (high) luminosity, the 
contribution of merging systems and/or AGN may be much higher. 
Moreover, at any specific frequency, the results here could be quite 
different - for example, in near and mid-IR wavelengths, AGN 
might be relatively much more luminous than cold dust emission 
from galaxy-wide starbursts (and even un-obscured AGN will con- 
tribute significantly at these wavelengths), thus the AGN luminos- 
ity density in such a rest-frame band might compete with or domi- 



nate t he luminosity density from star formation (see e.g. lBlain et al.1 
fl999T> . 

As discussed in § [2] there is an important difference between 
the total SFR density specifically induced via merger-driven galaxy 
starbursts (i.e. gas losing angular momentum owing to gravita- 
tional processes in the merger, falling to the galaxy center, and 
driving a short-lived starburst over ~ 10 8 yr timescales), and the 
total SFR density that might be identified observationally as "in 
ongoing mergers." The latter includes all star formation in systems 
that would be identified as merging (usually specifically limited to 
"major" mergers), via either some morphological or pair-separation 
based selection criteria. The duration of these phases (and hence 
the total SFR density associated with mergers in such a manner) 
depends on the exact selection criteria, but calibration of obser- 
vational method ologies with num erical simulations suggests it is 
~ 1 Gyr (see e.g. lLotz et alj |2008). During this time, except for the 
much shorter duration of the burst itself, the SFRs will (to lowest 
order, at least) reflect the "quiescent" or "normal" SFR of the disks, 
appropriate for their gas content and structu ral properties (see e.g. 
Idi Matteo etlu1l2007l ; ICox et alj|2008ll2009l) . In order to compare 
with these observations, we calculate the analogous SFR density 
in ongoing mergers with the following simple method: given the 
total rate of major mergers at each redshift, we simply assume a 
lGyr observable lifetime for each such (major, fi > 1/3) merger, 
and during this time assume it has a constant SFR equal to the rate 
of the quiescent systems with the same properties. We add the SFR 
density calculated in this fashion to that from the bursts themselves, 
and obtain an estimate of the "total" SFR that might be associated 
with e.g. disturbed or paired systems. 

This is much larger than the burst SFR density, especially at 
high redshifts. At high redshifts, merger rates are high, so a long 
observable duty cycle ~ 1 Gyr means that a large fraction of sys- 
tems will appear perturbed (i.e. the "merger fraction" will become 
large), and so a large fraction of star formation will appear in merg- 
ers. Here, we estimate this to rise from a few percent at z = to 
~ 20% at z ~ 2 and as high as ~ 20 — 50% at z ~ 3 — 6. Th i s is sim - 
ilar to the conclusions from the analysis in IHopkins et all J2006ch . 
using a different methodology but similarly attempting to calculate 
the total SFR in "ongoing" mergers. Of course, this should be larger 
than the burst SFR density at all times; but the primary reason the 
difference becomes so large at high redshift is that typical gas frac- 
tions are very large. As discussed in § [2] and shown in detail in 
simulations in lHopkins et alj ( l2009blli . large gas fractions lead to 
less efficient angular momentum loss and so (relatively) less effi- 
cient bursts in mergers, on average. However, local star formation 
in disks is not affected by this; star formation rates in the "quies- 
cent" or extended d isk mode continu e to rise super-linearly with 
/gas according to the lKennicuttl dl998l) law. 

3.6.2 Luminosity Densities: Comparison with Obsewations 

We compare these predictions to a number of observational con- 
straints. First, for the total luminosity/SFR density, we show 
the compilation of obser vations presented in IHopkins! d2004l) : 
IHopkins & Beacorrj J20061) . These come, for the SFR density, from 
a variety of different observations at various wavelengths; they are 
shown here in terms of the estimated total SFR density with the IR 
luminosity density following from the standard conversion adopted 
here. Complementary constraints at higher redshift can be inferred 
from observations of the Lym an-a forest in quasar spectra a nd ion- 
izing background, compiled in lFaucher-Giguere et alj j2008?). At all 
redshifts, the prediction is within the scatter of these observations; 
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Figure 8. Total IR luminosity density (and corresponding SFR density) as a function of redshift. We show the model prediction for the contribution from 
"normal" (non-merger-induced) disk star formation (green), and that for merger induced bursts (dark blue), and AGN (red), as in Figure|2] We also show the 
total IR luminosity density associated with "ongoing" interactions, which might be identified via morphology or pair-selected samples for a duration ~ 1 Gyr 
(light blue). Note that most of the star formation in these systems is the continuation of their "quie scent" star formation - the specifically merger-induced 
starburst lasts ~ 10 s yr. We compare the observational compilation from lHop kins & Bea coml l200d) for the total IR l umino sity density/SFR density (green 
diamonds) and at high redshifts the SFR density inferred from Lyman-Cf forest measurements in Faucher-Giguere et al. (2008, green circles). We also compare 
the results from the multi-wavelength "bolometric" AGN luminosity function compilation from lHopkins et alj l2007cl ) for the IR luminosity density in AGN 
(red). Observational estimates of the SFR density specifically induced b y mergers (e . g. sub tracting some "baseline" SFR estim ate from a control sample of 
non-merging systems for identified mergers) are shown (dark blue) from ljogee et alj 120091 triangles) and lRobaina et al.l i2009l inverted triangle); e stimates 
of the total SFR density in ongoing/id entifiable (e.g. pa ir or morphologically selected) mergers are shown (light blue) from lMenanteau et alj l200d circles), 
iBrinchmann et alj J 19981 squares), and lBell et alj feOOil star). At all redshifts, AGN represent a small contribution to the total (FIR-dominated) IR luminosity 
density. Mevger-induced star formation is a similarly small ~ 5% of the IR luminosity density; the total luminosity density associated with mergers is somewhat 
larger, but still small, rising from ~ 5 - 10% at z < 1 to ~ 10 - 20% at z ~ 2 - 4. 



at low redshifts z < 1, the median predicted is somewhat higher 
than the median of the observed points but the difference is small 
in absolute terms, ~ 0.2dex - well within the systematic uncertain- 
ties of both theory and observations. 

Next, we conside r the luminosity density in obscured AGN. 
iHopkins et al. (2007c) present bolometric quasar luminosity func- 
tions, compiled from observations at a wide range of different 
wavelengths, together with observationally inferred column den- 
sity distributions and template spectra. Adopting the fits therein, 7 
assuming that the obscured luminosity is re-radiated in the FIR, 
we construct the corresponding QSO IR luminosity density. This 
agrees well with our theoretical estimate. Similar constraints are 



obtained from the co mplimentary QLF compilations presented in 
IShankar et al.l d2009|). and from synthesis models of the IR back- 
grounds (e.g. lBlain et al.ll 1 999h . 

We also compare observational estimates of the luminos- 
ity/SFR density in mergers. First, several authors have attempted 
to estimate the total amount of star formation in observation- 
ally identified ongoing mergers or recent (morphologically dis- 
turbe d) merger remnants. We compare observation s compiled 
fromlBrinchmann et all i ll 9981) JlVlenanteau et alj j2006h : lBell et all 
fcooaT who estimate this quantity in morphologically-selected 
samples at z ~ — 1 .5. 8 Second, more recently, attempts have been 



7 A code for generating the observed quasar luminosity func- 
tions in various bands, based on these observations, is provided at 

http : / /www . cfa. harvard. edu/~phop kins /Site/qlf. html 



8 Note that most of these authors actually measure the fraction of the 
SFR density in or induced by mergers, not the absolute value. We con- 
vert this to an absolute density by rescaling with the observed total SFR 
density at the same redshift from the best-fit observed trend presented in 
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made to specifically isolate the merger-induced star formation rate 
density. Typically, in these cases, the SFR density of some merger 
sample (identified in a similar manner) is considered, but only af- 
ter subtracting away/removing the contribution from the expected 
"normal" mode star formation. In general, this is accomplished via 
comparison to some control sample of star forming galax ies with 
similar stellar masses and redshifts. iRobaina et al.l j2009h attemp t 
this from a pair-selected sample at z. ~ 0.4 — 0.8; |jogee et al.ld2009» 
consider a similar estimate from morphologically selected samples 
at z ~ 0.4— 1. Clearly, the two estimates (total SFR in "ongoing" 
mergers versus the SFR density enhancement from mergers) should 
be compared to the appropriate respective theoretical predictions, 
as discussed above. In both cases we see good agreement. At z ~ 1, 
the observations may in fact indicate the predicted growing differ- 
ence between "all" star formation in mergers and the star formation 
specifically induced by mergers. 



3.6.3 Contributions from LIRGs and ULIRGs 

In Figure [9] we illustrate the contributions to the luminosity den- 
sity from galaxies in various luminosity intervals - specifically, 
non-LIRG, LIRG, and ULIRG systems. For clarity we show just 
the results from our best-fit luminosity functions in Table [TJ the 
full allowed range scatters about these curves by ~ 0.15dex. We 
find the we ll-known result from a number of observational stud- 
ies (see e.g. lLe Floc'h et al 1l2005h : higher-luminosity systems pro- 
gressively dominate more of the IR luminosity density at higher 
redshifts. At z = 0, the luminosity density is dominated by rel- 
atively low-luminosity Lir ~ 10 10 Lq systems. The contribution 
from LIRGs (Lir > 10" Lq systems) rises rapidly from z = — 1, 
such that by z > 0.7 or so, these systems dominate the IR lumi- 
nosity density. Their fractional contribution remains relatively con- 
stant, above this point. ULIRGs (Lir > 10 12 Lq) also rise rapidly 
in prominence, from negligible contributions to the IR luminosity 
density (< 1%) at low redshifts to comparable ~ 20 — 50% contri- 
butions to LIRGs at z ~ 2, and by z > 3 dominating the total IR 
emission. The contribution from HyLIRGs (Lir > 10 13 Lq) rises in 
similar fashion, but is always much less than the contribution from 
ULIRGs - the highest contributions we find from such systems are 
at z > 3 at ~ 1 — 5% of the luminosity density. 

These trends simply reflect the predicted (and observed) evo- 
lution in the IR luminosity function break ~ L* (see Figure[3]l - i.e. 
the fact that all active systems become more IR luminous at high 
redshifts. It should be clear both from our discussion in § 13.31 and 
from direct comparison with Figure[8] that these luminosity classes 
do not necessarily represent distinct physical object classes. As we 
have shown, we expect most of the LIRGs at z ~ 1 , and ULIRGs at 
z > 2, to be "normal" star forming systems, with their luminosities 
driven by increasing gas content and rapid growth. 



4 CONCLUSIONS 

We present a simple model for the distribution of star formation 
rates and infrared luminosities owing to "normal" star-forming 
disks, merger-induced starbursts, and AGN. Comparing this with 
observations, we find reasonable agreement at z ~ — 3. At all red- 
shifts, we find that the low-luminosity population is dominated by 
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iHopkins & Beacom ( 2006j). Since this agrees well with our predicted total 
SFR density, it makes little difference if we use that instead. 



Figure 9. Contribution to the IR luminosity density from galaxies in differ- 
ent luminosity intervals. Top: Luminosity density, as Figure [8] We com- 
pare the total, and contribution from sub-LIRG (Lir < 10 11 Lq), LIRG 
(Lir = 10" - 1O 12 L ), and ULIRG (L m > 1O 12 L ) systems. HyLIRGs 
(Lir > 1O 12 L0) are negligible (< 1 — 5% contribution) in this total at all 
redshifts. For clarity, we show the results for the fits in Table[T] each curve 
has an approximate 0. 15 dex uncertainty. Bottom: Same, as a fraction of the 
total piR. LIRGs rise to dominance by z ~ 1. ULIRGs are comparable and 
then dominant in output at z ~ 2 and z ~ 3, respectively. We stress that, as in 
Figure|4] these luminosity cuts do not necessarily correspond to physically 
different classes of systems. 



disks, whereas the high-luminosity population becomes progres- 
sively more dominated by merger-induced bursts and then, ulti- 
mately, obscured AGN. 

The threshold for this transition is always at high luminosi- 
ties and low space densities. At higher redshifts, gas fractions in 
all systems increase - hence, specific star formation rates even in 
"typical" systems are much higher at high-z. As a consequence, the 
luminosity threshold for the disk-merger transition increases with 
redshift, from between LIRGs and ULIRGs in the local Universe, to 
ULIRG luminosities at z ~ 1 , to HyLIRG luminosities at z ~ 2 — 4. 
Similarly, the threshold for AGN dominance, at the bright ULIRG 
range at low redshifts, rises to HyLIRG luminosities at z ~ 1—2 
and to the very brightest HyLIRGs at z > 2. 

We provide simple fitting functions for each of these quan- 
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tities, and show how they depend on different parameters in the 
model. Most critically, it is the evolution in galaxy gas frac- 
tions that drives most of the evolution in the IR LFs. Obser- 
vations have shown that typical gas fractions in mas sive, star- 
forming galaxies increase rapidly with redshift (see e.g. Erb et al.l 
20061: iBouche et al.ll2007l;|Puech et al.ll2008l ; lMannucci et alj20(M 
Forster Schreiber et alj2009h . This naturally follows from the facts 
that cooling rates onto galaxies at high redshift are much higher 
than at low redshift, and there has simply been less time to pro- 
cess gas into stars (higher cosmic densities may also make stel- 
lar and AGN feedback relatively less efficient). A higher gas frac- 
tion by a factor of ~ 3, as implied by observations of Milky- Way 
mass disks at z ~ 2 — 3, leads to a factor ~ 5 — 6 higher star for- 
mation rate, according to the Kennicutt-Schmidt relation. At high 
masses/luminosities, i.e. where the number density of systems is 
falling exponentially, such a systematic increase in luminosity more 
than offsets the declining space density of massive galaxies with 
redshift. 

Other parameters have little effect. Changes in galaxy sizes 
and/or structural parameters, while potentially very important for 
e.g. how gas fractions are maintained, make little difference to the 
SFR distributions given some gas fraction and stellar mass distr ibu- 
tions. Allowing for a more steep index in the lKennicuttl d 19981) re- 
lation may help to account for the most luminous observed systems 
such as bright sub-millimeter galaxi es at the Hyper-LIRG threshold 
at z ~ 2 — 4 l lChapman et al.|[2005t) . However, the number counts 
of such objects remains quite uncertain, and recently it has been 
suggested that the average c ounts might be much lo wer with cos- 
mic variance still a concern dAustermann et al.ll2009l) . Larger sam- 
ples and better calibration of bolometric corrections - in particular, 
real knowledge of the appropriate dust temperatures for conversion 
to total-IR luminosities , which requires samp ling both sides of the 
cold dust peak (see e.g. I Younger et al ] |2009bl) - will be needed for 
better understanding of extreme systems. 

Our simple model succeeds reasonably well at explaining 
the observed global SFR density. At all redshifts, normal systems 
dominate the global SFR density. Obscured AGN contribute little, 
< 5%, to the total IR luminosity density. Substantial bias to IR- 
based SFR density estimates from obscured AGN would require an 
undiscovered population of heavily, isotropically obscured sources 
with luminosity densities ~ 5 times what is currently suggested 
(which would be in conflict with relic BH mass densities). 

The contribution of merger-induced bursts is similarly small, 
~ 5 — 10% at most redshifts. This owes both to the physics dis- 
cussed above - disks also rapidly increase their SFR density with 
redshift owing to higher gas fractions - and also to the fact that in- 
creasing disk gas fractions arbitrarily will not continue to increase 
the merger-induce d burst contrib u tion ar bitrarily. Rather, as dis- 
cussed in detail in iHopkins et al. I J2009bll3) . at high gas fractions 
angular momentum loss in mergers becomes less efficient - thus 
for an otherwise identical merger with a much larger gas fraction, 
the fraction of gas funneled into the nuclear starburst, relative to 
the total available, will be less (and the fraction that remains in 
an extended disk distribution to continue "normal" mode star for- 
mation will be larger), even if the absolute mass in the burst is 
larger. Similarly, at a given mass, the distribution of SFRs at all 
redshifts is dominated by normal-mode star formation - merger- 
induced bursts are important only in the high-SFR tail (~ 2 — 3 a) 
of the distribution. These trends explain a number of recent obser- 
vations that similarly indicate a small effect of merger enhance- 
ments to star formation rates, and that show that most star forma- 
tion by number and luminosity density appears to follow a simple 



trend or "main seque nce" as a function o f galaxy stellar mass and 
redshift felain et al.l [l999l: iNoeske et al Il2007bllil IPapovich et al.l 
l2006tlBell et alj|2005l : |jogee et al.ll2009l : lRobaina et alj2009l) . 

Support for these fractions also comes from completely in- 
dependent sources. Recently, a number of high-resolution stud- 
ies of spheroid formation via galaxy mergers have shown that 
properties such as the surface brightness profiles, sizes, concen- 
trations, kinematics, and isophotal shapes of spheroids are very 
sensitive to the mass fractions formed in such bursts, which pro- 
duce dense, disky, nuclear mass concentrations, versus the mass 
in a more extended envelope forme d via the violen t relaxation 
of the pre-burst stellar disks (see e.g. ICox et alj|2006l: iNaab et al .1 



| 200dlRobertson et alj2006l;lBurkert et alj2008l ; |jesseit et al.l2009l : 
IHopkins et al. 2008c a, l2009ch . In particular, typical ~ L„ early- 
type galaxies, which dominate the spheroid stellar mass den- 
sity, have properties that are reproduced accurately by simula- 
tions if and only if t his burst fraction is ~ 10% (for details, see 
IHopkins et alj|2009al le]). Independent analysis of their stellar pop 



ulation properties leads to similar conclusions (McDermid et al 



2009) 



20061 : ISanchez-Blazquez et alj|2007l : iReda et al .1120071 ; iFoster et al 



There is an important technical distinction between the total 
SFR density in "ongoing" or recent mergers and that actually in- 
duced by the merger (we present predictions for both). The former 
includes systems in their "normal" star-forming mode, observable 
for ~Gyr as perturbed or in pairs; the latter reflects specifically the 
~ 10 s yr event where gravitational torques drive a nuclear starburst. 
Under some circumstances, especially in very gas-rich mergers, the 
sum of this "normal mode" star formation over a long ~Gyr dura- 
tion yields significantly more total stellar mass formed than in the 
burst itself. The "ongoing" merger SFR density must, of course, rise 
with the observed merger fraction, reaching ~ 20% at z = 2 and as 
high as ~ 20 — 50% at z > 4; however we stress that these high 
fractions reflect predominantly the "normal" modes of star forma- 
tion simply present in systems that may be on their way to merging. 

Finally, we caution that the above comparisons are approx- 
imate, and intended as a broad means of comparing the primary 
drivers of star formation and their contributions relative to observed 
IR luminosity functions and SFR distributions. We have ignored a 
number of potentially important effects: for example, obscuration is 
a strong function of time in a merger, and may affect various lumi- 
nosities and morphological stages differently. Moreover, our simple 
linear addition of the star formation contribution of mergers to the 
IR LF and the AGN contribution is only technically correct if one or 
the other dominates the IR luminosity at a given time in the merger; 
however, there are clearly times during the final merger stages when 
the contributions are comparable. Resolving these issues requires 
detailed, time-dependent radiative transfer solutions through high- 
resolution simulations that properly sample the merger and quies- 
cent galaxy parameter space at each redshift, and is outside the 
scope of this work (altho ugh an import a nt subject for future, mo re 
detailed study; see, e.g. iLi et alj 2008; N arayanan et alj|2009b l l3; 
lYoungeretai]|2009al) . It would be a mistake, therefore, to read too 
much into e.g. the detailed predictions for sub-millimeter galaxies 
or other extreme populations in Figure [2] that may have complex 
dust geometries and/or a non-trivial mix of contributions from all 
of "normal" and "burst" mode star formation as well as AGN. How- 
ever, most of our predicted qualitative trends, including the evolu- 
tion of the luminosity density (and approximate relative contribu- 
tion of mergers) and the shift in where quiescent or merger-driven 
populations dominate the bright IR LF, should be robust. 
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Table 1. Fits to Model IR LF Predictions 



Object Class " 


L (±AL ) 6 


L' c 


L'"> 


4>o e </>'■' <t>" s 


«o" 


a" 


a"J 


ft* 


P" 


/?"<» 


Double Power-Law Fit (Full) - Standard Model 


Normal/Star-Forming 


11.37(0.14) 


2.17(0.82) 


-1.15(1.13) 


-3.97(0.50) 5.27(2.81) -6.11(3.53) 


1.57(0.45) 


-3.40(2.46) 


2.09(3.06) 


6.60(1.33) 


7.46(8.12) 


-17.57(10.81) 


Merger/Burst 


12.16(0.26) 


1.93(1.53) 


-2.41(2.09) 


-6.35(0.52) 5.15(3.06) -4.77(4.02) 


1.13(0.21) 


-0.86(1.30) 


0.30(1.76) 


3.23(0.59) 


3.94(3.55) 


-7.42(4.69) 


Obscured AGN 


12.15(0.29) 


2.07(1.63) 


-3.19(2.10) 


-6.62(0.47) 5.69(2.56) -5.39(3.15) 


0.73(0.23) 


-0.64(1.37) 


-0.83(1.85) 


2.55(0.50) 


3.90(2.91) 


-7.41(3.74) 










Double Power-Law Fit (First-Order in z) ■ 


■ Standard Model 










Normal/Star-Forming 


11.31(0.12) 


1.99(0.29) 





-3.30(0.41) 0.16(0.88) 


1.39(0.35) 


-1.68(0.72) 





6.59(1.10) 


0.95(3.44) 





Merger/Burst 


12.08(0.30) 


1.23(0.69) 





-5.67(0.59) 0.72(1.33) 


1.07(0.25) 


-0.53(0.55) 





3.18(0.60) 


0.79(1.71) 





Obscured AGN 


11.95(0.34) 


1.37(0.75) 





-5.80(0.53) 0.76(1.11) 


0.69(0.28) 


-0.83(0.60) 





2.40(0.45) 


1.06(1.32) 





Modified Schechter Function Fit - Standard Model 


Normal/Star-Forming 


10.46(0.58) 


2.35(0.97) 





-1.64(0.70) -1.28(1.14) 


0.46(0.74) 


-0.99(1.24) 





0.73(0.22) 


0.25(0.40) 





Merger/Burst 


10.23(1.83) 


1.05(0.64) 





-3.41(1.10) 0.21(0.57) 


0.28(0.72) 


-0.26(0.76) 





0.36(0.15) 


0.0(0.0) 





Obscured AGN 


10.10(0.66) 


0.90(0.60) 





-4.25(0.40) -0.36(0.81) 


0.0(0.0) 


-0.60(0.61) 





0.33(0.06) 


0.0(0.0) 






Double Power-Law Fit (Full) - Steeper Kennicutt-Schmidt Index 



Normal/Star-Forming 
Merger/Burst 
Obscured AGN 


11.28(0.14) 
12.39(0.25) 
12.15(0.29) 


2.43(0.84) 
2.06(1.49) 
2.07(1.63) 


-1.12(1.14) -4.09(0.53) 5.67(2.91) 
-2.59(2.07) -6.61(0.45) 5.00(2.65) 
-3.19(2.10) -6.62(0.47) 5.69(2.56) 


-6.35(3.58) 1.69(0.48) -3.93(2.55) 
-4.84(3.52) 0.96(0.17) -0.85(1.05) 
-5.39(3.15) 0.73(0.23) -0.64(1.37) 


2.48(3.09) 
0.48(1.43) 
-0.83(1.85) 


6.57(1.36) 
2.93(0.50) 
2.55(0.50) 


7.39(8.12) 
3.71(3.08) 
3.90(2.91) 


-17.75(10.66) 
-6.88(4.09) 
-7.41(3.74) 








Modified Schechter Function Fit 


- Steeper Kennicutt-Schmidt Index 











NormaVStar-Forming 10.34(0.56) 2.79(0.91) -1.62(0.69) -1.34(1.11) 0.51(0.71) -1.07(1.16) 0.72(0.21) 0.27(0.38) 

Merger/Burst 10.62(1.53) 1.05(0.59) -3.83(0.95) 0.30(0.51) 0.27(0.53) -0.21(0.59) 0.36(0.13) 0.0(0.0) 

Obscured AGN 10.10(0.66) 0.90(0.60) -4.25(0.40) -0.36(0.81) 0.0(0.0) -0.60(0.61) 0.33(0.06) 0.0(0.0) 



Double Power-Law Fit (Full) - No Gas Fraction Evolution 



Normal/Star-Forming 


11.38(0.15) 


1.06(0.96) 


-1.50(1.41) 


-4.00(0.57) 


4.36(3.45) 


-5.71(4.65) 


1.57(0.53) 


-2.12(3.36) 


1.05(4.67) 


6.70(1.52) 


7.34(9.66) 


-16.84(13.18) 


Merger/Burst 


12.20(0.26) 


1.35(1.54) 


-2.77(2.12) 


-6.40(0.53) 


4.51(3.18) 


-3.53(4.26) 


1.12(0.22) 


-0.21(1.42) 


-0.58(2.05) 


3.33(0.64) 


4.08(3.85) 


-7.61(5.03) 


Obscured AGN 


12.15(0.33) 


1.91(1.92) 


-2.88(2.62) 


-6.60(0.51) 


5.45(2.99) 


-5.29(3.90) 


0.75(0.24) 


-0.83(1.48) 


0.33(2.04) 


2.54(0.53) 


3.65(3.20) 


-6.66(4.23) 


Modified Schechter Function Fit - No Gas Fraction Evolution 


Normal/Star-Forming 


10.37(1.08) 


0.93(1.74) 





-1.52(0.94) 


-1.55(1.47) 





0.27(1.48) 


-1.01(2.53) 





0.69(0.34) 


0.18(0.61) 





Merger/Burst 


10.09(2.41) 


0.18(0.74) 





-3.41(1.12) 


0.58(0.61) 





0.18(0.92) 


-0.15(0.94) 





0.35(0.19) 


0.0(0.0) 





Obscured AGN 


10.07(0.71) 


0.73(0.68) 





-4.31(0.41) 


0.38(0.61) 





0.0(0.0) 


-0.28(0.65) 





0.32(0.06) 


0.0(0.0) 






s 

" Refers to sub-sample of objects for which the fit pertains. § 

* Parameters of best fit to the redshift-dependent form of the IR LF (Equations ll0ll2t . Lq is the break luminosity L* at z = 0, in log (Lq/Lq). f*j 

c ' d Redshift dependence of the break luminosity L* , per Equation lllK log {L* / Lq } = Lq + L' £ + L" £ 2 , where £ = log ( 1 + z) ■) ^ 

e Log LF normalization <f>* at z = 0, in Mpc~ 3 log - 1 Lir. 2 

f' s Dependence of normalization tj>* on redshift (log {(/>„ /Mpc~ 3 log - Lir} = </>o + <f>' £ + <j>" £ 2 ). 
h 

DenenHenre nf faint-end slone rv on redshift {rv — rvc -4- rv' P 4- r\ n £- 

l< 1 



Faint-end IR LF slope a at z = 0. r^J 

Dependence of faint-end slope a on redshift (a = ao + a' § + a" £ 2 ). c^' 
1 Bright-end IR LF slope j3 at z = 0. 

Dependence of bright-end slope f3 on redshift (/3 = /3 + j3' £ + fi" £ 2 ). §J 

O 

Observed galaxy mass functions extend to z ~ 4, the range used for these fits. Extrapolations beyond this redshift should be considered with caution. r^" 

If this is done, however, a minimum should be imposed on the bright-end slope at (i > 2. ^ 

to 
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APPENDIX A: PREDICTED LUMINOSITY FUNCTIONS 
AS A FUNCTION OF WAVELENGTH 

Figure [2] presents the predicted total IR luminosity functions from 
the models discussed here. To facilitate comparison with observa- 
tions (and broaden the range of observations to which we can com- 
pare), we here present corresponding predictions in a number of 
different rest-frame wavelengths. Figures IA1I IA2I IA3I IA4I & lA5l 
present the predicted LFs at rest-frame wavelengths of Sfi, 24 fi, 
60 /i, 100 /i, and 160 fi, respectively. For each, we compare to the 
available observations at or near that wavelength. 

For the star-forming systems (normal galaxies and mergers), 
we simply convert our predicted SFR and corresponding total IR 
luminosity to an observed luminosity in the given band, given 
the SED templates (themselves a function of bolometric luminos - 
ity) discussed in § 13.11 n amely those fro m Valiante et alj d2009l) . 
using the model SEDs in iDale & Heloul fe002h . As discussed in 
§ [3] varying the exact scaling of these corrections within observa- 
tional uncertainties is comparable to the uncertainty from adopt- 
ing different mass function estimators. For the AGN, we adopt 
the template SEDs f or obscured and unobscured systems from 
iHopkins et al. (2007c), but adopting alternative differe nt template 
obscured or unobscured AGN spectra (e.g. thos e in Elvis et al 
19941: IZakamska et alj |2004 IPolletta et alTl2006l ; iRichards et al 



warmer dust temperatures typical in AGN lead to their being rela- 
tively less important. 



2006al) makes little difference. 



We stress that these predictions should be regarded with con- 
siderable caution. The models here (or, for that matter, in any fully 
cosmological model for disk/merger/ AGN systems as a function of 
redshift) do not predict full SEDs. Rather, the robust quantity is 
some more physical number such as the total star formation rate 
(or AGN luminosity and obscured fraction). This allows robust es- 
timatiions of total IR luminosity, but we are now using a specific, 
simple empirical conversion between bolometric luminosity and lu- 
minosity at a given wavelength. If structural properties of galaxies, 
spatial distributions of gas and star formation, dust properties (gas- 
to-dust ratios, dust spatial distributions and dumpiness), clumping 
factors, and AGN contributions evolve, then these conversions will 
be problematic and may introduce systematic errors. In fact, it is 
very likely that these parameters that govern the SED do, in fact, 
evolve with redshift, or are different in merging and non-merging 
systems (given the different spatial distributions of gas and dust), 
and/or are a function of the relative AGN/star formation balance in 
the galaxy. 

More detailed modeling, including full, self-consistent radia- 
tive transfer treatment of high-resolution hydrodynamic simula- 
tions of mergers and normal galaxies with all the effects above in- 
cluded, will be necessary to predict e.g. the distribution of dust tem- 
peratures and other quantities critical, especially, for comparison 
with the number counts at long wavelengths (e.g. sub-millimeter 
galaxies). These models will be presented in future work (in prepa- 
ration); therefore we do not construct such comparisons (or attempt 
to compile predicted number counts) here. 

Nevertheless, the comparisons in Figures IAUA5I are infor- 
mative, and useful for future comparisons with observations, pro- 
vided appropriate caution is used. In particular, this allows us to 
see how the relative contributions of AGN and star formation vary 
as a function of IR wavelength (as the dust temperatures and SEDs 
are not the same). At shorter wavelengths, e.g. 8/z, AGN play a 
role at even moderate luminosities. We can compare to som e ob- 
servational studies, for example that in lBabbedge et alJd2006h . that 
explicitly separate the AGN and star-forming populations, and find 
good agreement. On the other hand, at the longest wavelengths, the 
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Figure Al. As Figure[2] but for the luminosity function in a specific rest-frame band (here, A = 8 (j.) as a function of redshift. We stress that we do not model 
the SEDs a priori, but simply adopt a specific set of empirical templates - as such, the information in this plot is identical to that in Figure [2] We compare 
to the same observations as Figure[2] for the observations at or near this rest-frame wavelength. Solid red pentagons show the estimates from lBabbedge et alj 
specifically for the contribution of AGN at this wavelength. 
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Figure A2. As Figure IaTI but at rest-frame 24/^. We compare observations spanning rest-frame 15 — 35 n (corrected with the same standard bolometric 
corrections to 24 n). 
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Figure A3. As Figure I A 1 1 but at rest-frame 60 fi. We compare observations spanning rest-frame 60 — 70 /j, (corrected with the same standard bolometric 
corrections to 60 /i). 
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Figure A4. As Figure IaTI but at rest-frame 100/.t. We compare observations spanning rest-frame 80— 120/i (corrected with the same standard bolometric 
corrections to 100/^). 
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Figure A5. As Figure lATI but at rest-frame 160 fi. We compare observations spanning rest-frame 120 — 350 fi (corrected with the same standard bolometric 
corrections to 160 (j,). 



© 0000 RAS, MNRAS 000, 000-000 



28 Hopkins et al. 



APPENDIX B: CONSEQUENCES OF MODEL 
ASSUMPTIONS 

The consequences of adding, removing, or changing some of our 
model assumptions are discussed in § 13.51 Here, in Figures lB"TllB2l 
IB3I & IB4I we explicitly illustrate the effects discussed there. 

In Figure [Bj] we reproduce Figure [2] but do not allow galaxy 
gas fractions to evolve with redshift (adopting the z — value at 
all redshifts). As discussed in § 13.51 this leads to significant under- 
prediction of the IR LF at high redshifts. In Figure lB2l we repro- 
duce Figure [2] again, but this time do not allow for disk sizes to 
be more compact at high redshift. This has a much smaller effect. 
In Figure IB3I we repeat Figure [2] but do not allow for scatter in 
any quantities (e.g. disk sizes, gas fractions, burst masses, quasar 
bolometric corrections, and SFRs at otherwise fixed properties); i.e. 
all values exactly trace the medians given in § [2] This suppresses 
the bright end of the LF. In Figure |B4 we repea t Figu re [2] but 



adopt a steeper power-law index for the iKennicuttl d 19981) relation 
(E» cx Egas), as discussed in § 13. 5t . These LFs correspond to the fits 
for the steep Kennicutt-Schmidt slope case presented in TableQ] 




Figure Bl. As Figure[2] but neglecting the increase in galaxy gas fractions 
with redshift. The SFRs of even normal, low-mass undisturbed disks are 
significantly under-predicted. 
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Figure B2. As Figure|2] but neglecting the evolution in disk sizes with red- 
shift. Because the observed size evolution (of star-forming galaxies) is rela- 
tive weak (R e (M t ) oc (1 + z)-(°~ 6 >) and the size evolution (at otherwise 
fixed properties) only enters into the SFR at sub-linear order, the difference 
is relatively small (factor ~ 2). 
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Figure B3. As Figure [2] but not allowing for any scatter in SFR at 
fixed galaxy mass, in mergers or disks (and no scatter in obscured frac- 
tions/bolometric corrections in AGN). The high-L tail is significantly sup- 
pressed. Note that the "kinks" at high redshift are artifacts of the analytic 
fitting functions used. 
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Figure B4. As Figure|2] but with all calculations adoptng a steeper index in 
the Kennicutt relation (Equation|3j, n s = 1.6, normalized to the same SFR 
for Milky Way-like disks. Systems at high redshifts are boosted significantly 
in SFR, and mergers are more concentrated in time, leading to sharper peak 
SFRs. The agreement with observations is somewhat improved, but overall 
the differences are comparable to the uncertainties from the adopted SFR. 
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